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This  study  presents  two  methods  of  analysis  for  the 
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including  resonance  and  a  unique  inflow  condition  are  also 
treated.   In  addition  a  nonlinear  method  of  characteristics 
solution  for  finite  cascades  is  described  which  permits 
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1.   INTRODUCTION 

1.1   General.   Unsteady  flow  in  turbomachines  can  be  broadly 
divided  into  two  categories :   circumferential  asymmetry  in 
the  flow  through  the  turbomachine  and  flow  perturbation  due 
to  the  self-excited  oscillatory  motion  (flutter)  of  the 
internal  blading  [1] .   As  shown  in  figure  1.1.1,  taken  from 
[2] ,  at  least  four  types  of  flutter  may  occur  in  a  compressor 
rotor.   The  two  stalled-f lutter  phenomena  are  encountered 
when  compressor  or  fan  blading  is  subjected  to  incidence 
angles  near  or  exceeding  the  stall  angle,  and  so  are  off- 
design  phenomena.   Subsonic  stalled  flutter  has  long  plagued 
the  engine  designer  [3] .   Choke  flutter  and  supersonic  unstalled 
flutter  are  low-incidence  instabilities  occurring  when  the 
flow  over  the  blade  is,  respectively,  transonic  along  most 
of  the  chord,  or  supersonic  over  the  outer  span  portions  [3]. 
Supersonic  unstalled  flutter  (hereafter  referred  to  simply 
as  supersonic  flutter) ,  which  has  appeared  only  recently  as 
a  significant  problem  in  turbomachinery ,  is  the  subject  of 
this  investigation.   Primarily  a  concern  in  compressors  and 
fans  [1] ,  it  is  an  instability  born  of  the  recent  technologi- 
cal advances  aimed  at  increasing  the  thrust-to-weight  ratios 
of  turbomachines:   increased  tip  speeds,  composite  blades, 
lighter  discs  and  blades,  and  the  elimination  of  part-span 
shrouds  [2,3].   Some  insight  into  the  nature  of  supersonic 
flutter  is  gained  by  a  consideration  of  its  characteristics: 
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Figure  1.1.1   Compressor  map  showing  boundaries  for  four 

types  of  flutter. 


1.  Its  flutter  boundary  imposes  a  high-speed  limit  on 
compressor  operation  [2]. 

2.  It  is  a  low-incidence,  attached-flow  instability  that 
can  occur  at  the  design  condition  (figure  1.1.1). 

3.  Flutter  amplitude  decreases  with  increasing  back 
pressure  on  the  rotor  [2] .   (For  this  reason  the  present 
analysis  considers  the  most  critical  case  of  zero  pressure 
rise  across  the  rotor.) 

4.  At  the  flutter  condition  all  the  blades  are  observed 
to  flutter  at  the  same  frequency  with  a  constant  interblade 
phase  angle  [2] . 


5.  Both  single-degree-of -freedom  torsional  flutter  (pitch) 
and  two-degree-of-freedom  torsional-f lexural  flutter  (pitch 
and  plunge)  can  occur  [3,4]. 

6.  The  flutter  boundaries  are  significantly  affected  by 

airfoil  shape  [4,5]. 

7.  The  primary  flutter  mechanism  is  an  inviscid  phenomenon 
If  the  unsteady  aerodynamic  force  lags  the  blade  motion 

by  a  phase  angle  which  is  sufficient  to  ensure  positive 
work  over  the  entire  cycle,  then  flutter  can  occur  [3,4]. 

It  is  clear  that  reliable  techniques  for  predicting  the 
unsteady  pressure  distribution  on  a  fluttering  blade  are 
required  if  one  hopes  to  predict  and  avoid  the  onset  of 
flutter.   With  this  end  in  mind  the  cascade  has  been  exten- 
sively employed  as  a  theoretical  abstraction  of  a  compressor 
rotor.   The  unsteady  flow  is  considered  in  an  annulus  of 
differential  radial  height  within  which  equispaced  airfoils 
are  radially  disposed.   Unwrapping  such  an  annulus  gives 
rise  to  the  two-dimensional  approximation  of  plane  flow 
through  a  cascade  of  airfoils  [1] .   The  cascade  is  termed 
linear  if  the  airfoils  are  flat  plates  of  zero  thickness. 
As  seen  in  figure  1.1.2,  a  rotor  does  not  require  a  super- 
sonic axial  flow  to  have  a  supersonic  velocity  component 
relative  to,  and  parallel  to,  the  blade  axis.   Indeed  this 


XR  A  Arnold!  et  al.  in  "Supersonic  Chordwise  Bending 
Flutter  in  CascadesT-^ratt  *  Whitney  Repor  o  •  ^-5271 
of  31  May  1975  show  that  chordwise  bending  can  also  occur. 
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Figure    1.1.2      Supersonic   cascade   classification 


distinction  provides  a  convenient  method  of  classifying 
supersonic  cascades.   If  the  axial  velocity  is  subsonic, 
the  cascade  is  said  to  have  a  subsonic  leading  edge;  in  the 
other  case  it  has  a  supersonic  leading  edge  [6] .   The  physi- 
cal significance  of  this  distinction  is  clear  from  figure 
1.1.2.   In  the  supersonic-leading-edge  cascade  each  blade 
influences  only  its  immediate  neighbors.   When  the  cascade 
has  a  subsonic  leading  edge,  however,  disturbances  generated 
by  the  upstream  portions  of  each  blade  affect  the  flow  around 
each  subsequent  blade  in  the  cascade;  furthermore  the  flow- 
field  at  the  underside  of  each  blade  is  influenced  by  the 
wake  disturbances  of  all  preceding  blades.   Since  the  supersonic- 
leading-edge  cascade  has  been  thoroughly  investigated  [7,8,9, 
10],  the  more'  difficult  subsonic-leading-edge  case,  which  is 
also  the  more  interesting  in  practice,  is  considered  in  the 
present  study.   Another  cascade  classification  is  frequently 
employed:   if  each  blade  is  assumed  to  be  preceded  and  followed 
by  infinitely  many  others,  the  cascade  is  said  to  be  infinite; 
if  there  is  an  identifiable  first  blade  —  a  necessary  pre- 
requisite for  many  numerical  solutions  —  the  cascade  is 
termed  finite. 

A  considerable  amount  of  effort  has  been  expended  in 
attempting  to  predict  —  numerically  and  analytically  — 
the  unsteady  pressure  distribution  over  a  blade  in  an  oscilla- 
ting supersonic  cascade  with  subsonic  leading  edge.   A 
brief  discussion  of  some  of  the  more  noteworthy  advances 
follows . 


1 . 2   The  Linear  Subsonic-Leadinq-Edge  Cascade.   Verdon 
and  McCune  [11]  formulated  a  linear  boundary  value  problem 
for  the  subsonic-leading-edge  configuration,  solving  it 
for  velocity  potential.   This  permitted  calculation  of  the 
pressure  along  the  entire  upper  surface  of  an  arbitrary  blade 
and  along  the  lower  surface  from  the  leading  edge  to  the 
point  of  impingement  of  the  trailing-edge  Mach  wave  from 
the  preceding  blade.   In  a  more  recent  paper  Verdon  [12] 
formulated  a  second  boundary  value  problem,  in  which  pressure 
is  the  basic  dependent  variable,  which  permitted  determination 
of  the  pressure  along  the  entire  lower  surface  as  well. 
Nagashima  and  Whitehead  [13]  and  Goldstein  [14]  obtained 
solutions  by  replacing  the  cascade  blades  with  doublet  dis- 
tributions.  Kurosaka  [6]  introduced  what  he  terms  the 
passage  approach  in  applying  Laplace  transform  techniques 
to  a  slowly-oscillating  cascade,  obtaining  closed-form 
analytical  expressions  for  the  pressure  distributions  of 
certain  configurations.   He  recently  [15]  extended  his  solu- 
tion to  include  arbitrary  frequencies  of  oscillation.   In 
assuming  that  a  space  periodicity  of  the  flow  exists  from 
one  blade  passage  to  the  next,  Kurosaka  claims  that  he 
avoids  complications  arising  from  the  breakdown  of  linear 
acoustic  theory  in  the  far  field  [4] .   Platzer  and  Brix 
[16]  and  Platzer  and  Bell  [17]  employed  the  method  of  char- 
acteristics to  solve  directly  for  the  perturbation  veloci- 
ties and  pressures  everywhere  in  the  flowfield.   Their 
computer  program  has  been  thoroughly  checked  against  existing 


Figure  1.2.1   Limits  of  applicability  of  Chadwick's  solution 


solutions  [11,12,13]  and  is  extensively  employed  in  the 
present  study  as  a  method  of  comparison  for  both  the  analy- 
tical and  numerical  solutions  presented.   Chadwick  [5]  used 
Sauer's  solution  [18]  of  the  unsteady  potential  equation, 
accurate  to  a  first  order  in  frequency  of  oscillation,  to 
determine  the  pressure  distribution  on  a  given  blade  from 
the  leading  edge  to  the  point  of  the  first  reflection  (point 
A  in  figure  1.2.1)  of  disturbances  from  the  preceding  blade. 
He  obtained  exact  analytical  agreement  with  Kurosaka  [6]  over 
that  portion  of  the  blade  for  the  single-degree-of-f reedom 
pitching  cascade.   Chadwick's  technique  has  the  decided 
advantage  of  being  much  simpler  than  that  of  Kurosaka,  easily 


yielding  closed-form  expressions  for  any  cascade,  and 
raising  the  question  as  to  whether  the  method  could  be 
extended  to  the  entire  blade  for  the  cases  of  both  pitch 
and  plunge.   Finally,  Adamczyk  [19]  ,  in  a  yet-unpublished 
work,  employs  the  Wiener-Hopf  method  to  obtain  the  pressure 
distribution  for  the  linear  cascade  with  subsonic  leading 
edge. 

1 . 3   The  Nonlinear  Subsonic-Leading-Edge  Cascade.   Attempts 
to  deal  analytically  or  numerically  with  an  oscillating 
cascade  of  blades  having  thickness  are  not  nearly  so  numerous 
as  their  linear  counterparts.   Carrier  [20]  provided  an  exact 
theory  for  a  single  oscillating  wedge.   VanDyke  [21,22]. 
also  .treated  the  isolated  wedge  as  well  as  an  oscillating 
airfoil  of  arbitrary  profile.   Teipel  [23,24]  employed  an 
ingenious  method-of-characteristics  approach  to  numerically 
obtain  the  unsteady  pressures  on  a  single  oscillating  profile 
of  arbitrary  shape,  obtaining  good  agreement  with  the  single- 
blade  solutions  of  Carrier  and  VanDyke.   Teipel' s  method 
assumes  potential  flow  with  negligible  entropy  changes 
across  leading-edge  and  trailing-edge  shocks.   Kurosaka  and 
Edelfelt  [4]  have  obtained  numerical  results  for  an  isolated, 
oscillating,  parabolic-arc  airfoil  and  are  presently  working 
to  extend  their  solution  to  a  cascade  of  such  airfoils. 
Chadwick  [5]  employed  Teipel 's  method  to  obtain  the  unsteady 
pressure  distribution  over  the  preinterf erence  zone  of  an 
arbitrary  blade  in  a  cascade  of  wedges.   His  results  indicate 


that  nonlinear  thickness  effects  significantly  alter  the 
unsteady  pressure  profile  of  an  arbitrary  blade. 

1.4   The  Present  Investigation.   The  present  work  has  two 
major  parts:   The  first  is  analytical  and  treats  the  linear 
cascade  while  the  second  is  numerical,  dealing  with  a  non- 
linear configuration. 

In  the  linear  effort  a  closed-form  analytical  solution 
for  the  pressure-difference  distribution  along  a  blade  in 
a  slowly-oscillating  cascade  is  obtained.   The  solution 
technique  is  that  developed  by  Chadwick  [5]  in  the  pre- 
interference  zone  (figure  1.2.1)  of  a  pitching  cascade,  in 
which  Sauer's  solution  of  the  unsteady  linearized  potential 
equation  is  employed.   In  the  present  work  Chadwick 's  solu- 
tion is  extended  through  all  reflection  zones  and  into  the 
blade's  wake.   In  addition,  the  solution  is  obtained  for 
two  degrees  of  freedom  (pitch  and  plunge) .   Two  standard 
configurations  —  Verdon's  cascades  A  and  B  [12]  —  are  con- 
sidered, one  from  a  finite-cascade  approach  and  the  other 
as  an  infinite  cascade  in  which  Kurosaka's  passage  approach 
is  employed.   Extensive  comparisons  are  made  with  existing 
solutions  [6,12,19,17];  and  several  special  topics,  including 
cascade  resonance  and  a  unique  inflow  condition,  are  briefly 
treated. 

The  nonlinear  effort  is  an  evaluation  of  the  feasibility 
of  extending  Teipel's  unsteady  method  of  characteristics 
to  a  cascade  of  airfoils  with  arbitrary  profiles.   Teipel 


[23,24]  represents  the  nonlinear  flowfield  by  a  harmonic, 
unsteady,  velocity  potential  superimposed  on  a  steady  one, 
and  discontinuities  in  the  flowfield  variables  across  the 
Shockwaves  emanating  from  the  airfoil  are  determined  from 
the  isentropic  Rankine-Hugoniot  relations.   These  concepts, 
coupled  with  the  characteristic  grid  and  compatibility  relations 
derived  by  Teipel,  permit  determination  of  the  steady  and 
unsteady  preinterference  pressures  acting  on  the  second 
blade  of  a  cascade  of  parabolic-arc  airfoils.   The  numerical 
results  are  compared  with  experimental  data  from  tests  con- 
ducted with  a  pitching  cascade  by  Fleeter  and  Riff el  [25J 
and  with  linear  theory.   An  improved  method  for  determining 
the  jump  conditions  across  the  bow  shock  is  presented  and 
'the  computational  procedure  is  carefully  outlined  to  facili- 
tate extension  of  the  technique  to  the  entire  flowfield  of 
a  completely  general  cascade. 

1.5   Some  Comments  on  Notation.   In  the  linear  portion  of 
this  investigation  it  is  the  unsteady,  dimensionless , 
pressure-difference  distribution  along  a  particular  cascade 
blade  that  is  sought.   This  pressure  difference  is  a  complex 
quantity  given  the  symbol  AC   and  represents  an  amplitude 

ir 

function  of  distance  along  the  blade.   The  harmonic  oscilla- 

ikt 
tion  factor   e     is  always  implicitly  understood.   In 

addition  it  is  desirable  to  obtain  an  expression  for 

AC    at  the  instant  that  the  blade  in  question  achieves 
p 


maximum  up  displacement,  referred  to  as  the  initial  position 

This  expression  is  obtained  by  multiplying  through  by  the 

factor   e   '   and  so  it  is  indicated  by  the  sumbol  Ac  e  / 

P 

In  the  nonlinear  portion  of  this  study  only  the  local 

pressures  —  C   upper  and  lower  —  are  sought.  They  are 

also  complex  quantities  and  it  is  convenient  to  express 

them  as  a  magnitude  and  phase  angle  relative  to  the  motion 
of  the  blade  under  study;  the  symbols 
respectively,  are  employed. 


C 

P 


and   /C 

P 
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2.   THE  LINEAR  LOW-FREQUENCY  CASCADE 

2.1  General.   A  cascade  of  flat  plates  with  chord   c 
oscillating  harmonically  with  frequency  G>   in  a  steady, 
uniform  flow   U   is  said  to  oscillate  at  a  dimensionless 
reduced  frequency  of  oscillation  k  defined  by 

*=-^J-  (2.1.1) 

2 

If   k   is  small  when  compared  with  k  ,  then  the  cascade  is 

termed  a  low-frequency  or  slowly-oscillating  cascade.   Although 
low-frequency  flutter  is  not  common,  it  is  worthy  of  investi- 
gation for  several  reasons,  the  major  one  being  that  it  is 
the  only  case  for  which  one  can  hope  to  attain  a  relatively 
simple,  closed-form  expression  for  the  unsteady  pressures 
on  a  given  blade.   Such  expressions  provide  insight  into  the 
flutter  phenomenon,  as  evidenced  by  the  identification  of 
a  unique  inflow  condition  in  a  subsequent  chapter  of  this 
study.   Also,  the  low-frequency  cascade  provides  an  important 
limiting  case  for  any  general-frequency  solution,  numerical 
or  analytical. 

2. 2  Problem  Formulation. 

2.2.1   The  Linearized  Potential  Equation.   The  flow 
variables  in  the  cascade  flowfield  are  assumed  to  be 
composed  of  a  freestream  component  and  a  much  smaller 


perturbation  component 

P=FU-vP'  (2.2.1) 

The  cascade  is  assumed  immersed  in  a  perfect  gas,  each 
particle  of  which  maintains  constant  entropy  as  it  moves 
through  the  flowfield 

•||-=0  (2.2.2) 

In  addition  the  flow  is  assumed  inviscid  and  irrotational, 
guaranteeing  the  existence  of  a  velocity  potential 

V=V#  (2.2.3) 

The  continuity  and  momentum  equations  for  the  flow  are, 
respectively 


& 


Dt- 


'  P   V- 


•  V  =  O         (2.2.4) 


DV     7P    — 

=  0  (2.2.5) 


introducing  the  velocity  potential    and   integrating    (2.2.5) 
produces,    after   linearization 

I  J^+XjjlHL     +     1^  =  0  (2.2.6) 

at.  ^x  />«> 

Using  the  equation  for  the  freestream  speed  of  sound 

C^^V^r-  (2.2.7) 

the  continuity  eouation  can  be  linearized  to  give 

3bL  +.  ^  +  -J r^L.  ^. -u- il-1  s0    (2.2.8) 

***       *<?1       A»  *i     *>*•  ** 
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From    (2.2.6)    one   obtains 

St 

and 


z?  .    r^tf       —  *Na 


9X 


pSSt^TSr]   (2'2-10 


) 


Substitution  of  (2.2.9)  and  (2.2.10)  into  (2.2.8)  then 
produces  the  linearized  unsteady  potential  equation  in 
two  dimensions 

2£w  ^i;S(il5  "'"  }75?  "**  T^r-    (2.2.11) 

where  M  is  the  free  stream  Mach  number 

M    =  -LL-  (2.2.12) 


2.2.2   The  Linearized  Pressure  Coefficient.   A  generalized 
form  of  Bernoulli's   eauation  follows  directly  from  the 
momentum  equation  (2.2.5) 
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lt"+   t   +  J"7S  CONST.       (2.2.13) 


Introducing  the  perturbation  quantities  in  the  form  of 
(2.2.1)  and  ignoring  the  products  of  perturbation 
quantities  one  obtains 


f^^-^O        (2.2.14, 


where  the  prime  notation  has  again  been  dropped.   Thus,  the 
dimensionless  pressure  coefficient  is  given  by 

2.2.3   The  Linearized  Boundary  Conditions.   The  surface 
of  a  body  placed  in  the  flow  can  be  described  as  a  function 
of  space  and  time 

n^UrO  (2.2.16) 

where  r   is  the  position  vector  of  points  on  the  surface. 
The  condition  that  the  flow  at  the  body  surface  be  everywhere 
tangent  to  the  surface  can  be  expressed  by 


-LEL  -o  (2.2.17) 

•pt. 
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The  plunging  flat  plate. 

The  pitching  flat  plate. 

Figure  2.2.1   Two-degree-of-freedom  oscillation. 

at  the  surface.   If  the  body  is  a  flat  plate  oscillating 
harmonically  in  plunge  and/or  pitch,  as  shown  in  figure  2.2.1, 
then  the  describing  functions  are  respectively 


Cu»"t 


F»<W=rh'e    s° 


(2.2.18) 


'5^^V  +  0''W"We   a°       (2.2.19) 


where   h    and  <*    are  the  maximum  displacements  and   b 
o        o 

is  the  pitch  axis  location  on  the   x   axis.   Applying  the 
generalized  flow  tangency  condition  (2.2.17)  to  (2.2.18) 
and  (2.2.19)  yields  the  boundary  conditions  for  plunge  and 


pitch   respecitvely 

* 

|^-  =    -hoicoeiwt      cL-t    ^  =  o  (2.2.20) 

-|£-  =   -ex^Tl-v  Lco(x-b)]ec,>>,t    o.t^:0        (2.2.21) 

2.2.4   Nondimensionalization.   It  is  convenient  to  intro- 
duce nondimensional  variables  into  the  governing  equations 
to  simplify  and  generalize  them.   The  distance  variables 
(x,y,b)   are  normalized  by  chordlength   c  ,  velocity  potential 
by  the  product   Uc   and  time  by  the  quantity   c/U  .   It  is 
also  convenient  to  set 

h„=  c      **:=  1  RADIAN        (2.2.2-2) 

and  to  separate  the  velocity  potential  into  a  complex 
amplitude  and  phase  angle 

tfcx^/fc^  9*(*/*>ei>"t        (2.2.23) 

Thus  the  potential  equation  (2.2.11),  pressure  coefficient 
(2.2.15)  and  boundary  conditions  (2.2.20,21)  become 

a1^.  -  £SL   -**Mz0-*-;u*M2!4  =  O    (2.2.24) 
Cp=^2(i^ i*$)  (2.2.25) 


—  =  -v.fe    a^  ^so        (2.2.26) 


■55-=  —  [!"»-  ^CX-bJJ  oit^o      (2.2.27) 


where , 


/ 


Z^MZ-1  (2.2.28) 


and   C    and  ^4  are  complex  amplitude  functions  of   x   and 

ikt 
y   only  with  the  factor   e     omitted. 


2.2.5   Cascade  Geometry.   The  physical  parameters  chosen 
to  describe  the  cascade  are  the  vertical  distance  between 
the  blades   d   and  the  horizontal  distance  between  the 
leading-edge  Mach  waves   m  .   Both  parameters  are  normalized 
by  chordlength.   The  cascade  can  also  be  described  by 
solidity   £   and  stagger  angle  ©  ,  which  are  shown  in 
figure  2.2.2.   These  parameters  are  related  as  follows 


T*N@=-*T^T  (2.2.29) 


j&»  -S^L  (2.2.30) 


where  &     is  given  by  (2.2.28).   The  linear  cascades  employed 
in  this  study  are  Verdon ' s  cascades   A   and   B   [12],  which 
have  become  standards  in  the  field  and  are  depicted  in 
figure  2.2.2. 


2.3   The  Isolated  Airfoil.   Consider  a  single  flat  plate 
oscillating  with  two  degrees  of  freedom  as  shown  in  figure 
2.2.1,  and  define 

S^*-^    s^-K-t-yjj;  (2.3.1) 

Then   s   is  constant  along  any  left-running  (upper  surface) 
Mach  wave,  while   s   is  constant  along  right-running  (lower 
surface)  waves.   Sauer  [18]  showed  that  the  potential 
equation  (2.2.24)  has  the  following  left-running  and 
right-running  solutions,  respectively,  for  a  slowly- 
oscillating  airfoil 


#^->^  +  *[N(^-  *jr-29*<$        (2.3.2) 

ftl%$=$u<3\+  *\jnjX*+  ^^  £ut^  (2.3.3) 


where   the      g      and      h      functions    are   determined   from  the 
boundary   conditions    (2.2.26,27).      Requiring  that   the 
velocity   potential   be   continuous   across   the   leading-edge 
Mach  waves,    one  has 

£u(o}=  ViuCc»  =S  £u<o>  -Mo)  -  O  (2.3.4) 

Enforcing  flow  tangency  at  the  upper  surface  of  the  airfoil 
yields,  for  plunge  and  pitch  respectively 


Cascade  parameters. 


M:   1.345 

p:   0.899 

m:    0.319 

d:    0.400 

£z    1.268    (1/D 

0:    30.5° 


pitch  axis 
at  midchord 


M:  1.281 

£z  0.801 

m:  0.360 

d:  0.301 

ii  1.490    (1/L) 

©:  26.6° 


pitch  axis 
at  midchord 


Figure    2.2.2      Cascade   geometry, 


PL:       ^uCS>  =  o       huC&>=-L.S         (2.3.5) 


PI:      S«{%^jr         Wul«  — ^(&.^byi"«)  (2.3.6) 


And  so  on  the  upper  surface 


PL:  #<*/*>  *^&  (2.3.7) 

PI:      ^■»/^«y-^[r^&Cb^yaM^    (2.3.8) 

Thus  the  upper-surface  pressure  coefficient  follows  from 
(2.2.25) 


PL:  -'      c£  =  -^£  (2.3.9) 

PI:  Cn»-7"+^[*l|-^}+V']     (2.3.10) 

where  the  subscript  indicates  the  single  airfoil.   Similarly, 
for  the  lower  surface 

C^s-Cjf,  (2.3.11) 

2. 4   Continuity  of  the  Velocity  Potential  Across  a  Mach  Wave. 
Consider  two  reflection  zones  separated  by  the  Mach  wave 
s  =  s   ,  as  shown  in  figure  2.4.1,  and  let  the  variables   z 
and   z   respectively  measure  distance  tangent  and  normal  to 
the  Mach  wave.   The  velocity  potential  in  zone  II  is  simply 


ZONE  I  /  ZONE  II 


Mach  wave:  s  =  s, 


z 


-*-x 


Flat-plate  airfoil 


igure  2.4.1   Variation  of  the  velocity  potential  along 
the  boundaries  of  reflection  zones. 


at  in  zone  I  plus  an  additional  disturbance 


$*.*»+  ^ 


(2.4.1) 


nee  the  tangential  velocity  component  remains  unchanged 
ross  the  wave,  (2.4.1)  implies  that  $d   is  a  function  of 
alone,  along  the  Mach  wave 

q^s^CZ)    4-oy   S^So      (2.4.2) 

it   since      z      is    itself   constant   along   the   Mach  wave,    then 
must  have   the   form 


<£,t=.  CONST.         -for     s=s>a 


(2.4.3) 


That  the  constant  is  zero  follows  from  the  physical  con- 
straint that  the  velocities  remain  finite  in  the  wave 
itself,  and  so 

^CSo^s  <^B  CS0>  (2.4.4) 

which  is  in  agreement  with  the  conclusion  reached  by 
Landahl  [26] . 

2 . 5   The  Finite  Two-Degree-of-Freedom  Cascade.   As  an 
elementary  means  of  introducing  the  method  of  analysis, 
the  pressure-difference  distribution  is  found  for  the  second 
blade  of  a  finite  version  of  cascade  A.   This  analysis  also 
affords  the  opportunity  to  compare  the  finite-cascade  .results 
obtained  with  numerical  finite-cascade  data  [17]  and  analy- 
tical results  [6]  of  existing  solutions.   It  is  emphasized 
that  the  pitch  results  through  equation  (2.5.27)  were  first 
obtained  by  Chadwick  [5].   The  reflection  zones  referred  to 
are  pictured  in  figure  2.5.1,  and  the  pressure  results  are 
given  in  terms  of  the  isolated-blade  pressure  coefficients 
(2.3.9,10).   Before  proceeding  with  the  analysis,  some 
additional  variables  are  introduced  for  blade  two 

tr=^-ya^       ^  =  ^y&"X      (2.5.1) 

which  are  related  to  those  of  blade  one  by  the  translation 
relations 


Cascade  A  (not  to  scale) 
Mach  number:  1.345 


Figure  2.5.1   The  finite-cascade  reflection  zones. 


x=  $  -»•  m-^a     ^=  i  -v-<J 


(2.5.2) 


s=  cr-»-  m    s  =  a- -*-  vn  -*•  a 


Lysa 


(2.5.3) 


2.5.1   Zones  I  and  IA.   In  zones  I  and  IA  Sauer's  single- 
blade  solution  applies 


ti!**tlfl   =  ^<&>^*[*«^-^c7c?»(s>]     (2.5.4) 
$?%V  =  ^a<^  +  *£^>-^>^>]   (2.5.5) 


and  in  zone  IA  the  boundary  conditions  take  the  form 


PL:         -^2^-  =-tfee1/*    <x.t    ^=0        (2.5.6) 


-  -M-i-LteCQr-bile^  <Tt  »^0   (2.5.7) 


where   m  is  the  interblade  phase  angle  by  which  the  second 
blade  leads  the  first.   Applying  the  boundary  condition  at 
the  upper  surface  of  each  blade  and  using  that 

£,(crt  =  Vi.to^  •=.   O  (2.5.8) 

gives,  for  plunge  and  pitch 

PL:  £V(S>=0       ^AC<T)=A        (2.5.9) 


PI: 


£,<«  =  -£-     ^fcn^^Le^nrB    (2.5.10) 


PL:      h.lSO^^rS     hv*(<r}=:i-<reV*+C!    (2.5.11) 

r  /* 


PI:  P  (2.5.12) 


where   A,  B,  C   and   D   are  constants  of  integration.   By 
enforcing  continuity  of  the  velocity  potential  across  the 
leading-edge  Mach  wave  of  blade  two,  the  constants  are 
evaluated  and  the  velocity  potentials  determined  in  terms 
of  the  isolated  blade  potential  <j>x   (2.3.7,8) 


PL:     ^-^f         0*\s#e>-*  it™     (2.5.13) 

r 


PI:  (2.5.14) 

where   £   is  a  constant  that  need  not  be  evaluated  because 

2 

it  becomes  the  coefficient  of  a  term  in   k   when  the 

pressure  coefficients  are  formed.    Thus,  on  the  upper 
surface  of  the  second  blade,  in  terms  of  the  isolated-blade 
pressure  coefficient   Cp(  (2.3.9,10) 

PL:  Cp2**=s  C£e'^  (2.5.15) 


PI:  c£x*^  C^e^  +  5^m  (2.5.16) 


2.5.2   Zones  IIB  and  VIB.   In  zone  IIB  the  velocity 
potential  is  that  of  the  isolated  blade  (2.3.3) 

0f8=  <j>\  (2.5.17) 

while  in  VIB  the  potential  acquires  an  additional 


"^The  symbol   £   is  used  throughout  section  2  to  indicate 
a  constant  which  need  not  be  determined. 
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right-running  disturbance 


0"!**=  &**   +^6CS>  +  te[h6(S^  ^xV^tS>]         (2.5.18) 

2.5.3  Zone  II.  In  zone  II  a  right-running  reflection 
occurs  and,  using  the  translation  relations  (2.5.2,3)  to 
shift  to  the  blade- two  coordinate  system,  one  has 

PL:      if  =  9^?^z(^>^Kfe  ^b2(^)^  y^(T-)]     (2.5.19) 


PI:  (2.5.20) 


-^*(5*)  +■  *[h2tcf)  ^  ^  t  £*  <«] 


where 


^*(o*  =t Viz(o^  « O  (2.5.21) 

Applying  the  flow  tangency  condition  at  the  lower  surface 
of  blade  two  and  using  the  Sauer  functions  to  construct  the 
velocity  potentials  yields 

PL:  $?  =  0f^^-fr(v-e;-/4)  +  Vnl       (2.5.22) 

p   L  J 

PI:        -v  ^[-^am  +>SciMl+  Cb^-^M^Kl-e'v*^     (2.5.23) 


And  the  pressure  coefficients  follow  from  the  potentials 


*-n 


PL:  cfz  -^(2-e^)  (2.5.24) 


PI:         Cpf  as.  Cp^Z-e^i-i-^Cm-Zfi^-^adM*)  (2.5.25) 


Thus  the  pressure  difference  can  be  computed  through  zone 
II,  and  by  multiplying  by  the  factor  e  /  the  pressure- 
difference  amplitude  is  obtained  at  the  moment  the  second 
blade  reaches  maximum  displacement 

PL:   ^Cpe-^^-zcfto-e-v*)     o±  ^   <  \-m  (2.5.26) 


PI:  (2.5.27) 


2.5.4   Zone  III.   In  zone  III  a  left-running  wake  distur- 
bance is  experienced,  and  after  applying  the  translation 
relations 


PL: 


fi*   -  *?+f*W+*[ts&-^J3i&  +  *]  (2.5.28) 


PI:       -s(y^^y6<lMa)^[-^^M^(2^\iMx+flV»iMl)Cl-e^)]  (2.5.29) 
*#1  ^JCS)  +Vf  [hj(s^-.ilp^^j(s>] 


where 


^jtO  -  In,  Cl)  =  Q  (2.5.30) 

Since  zones  III  and  VIB  are  separated  by  the  wake  of  the 
blade  rather  than  the  blade  itself,  the  wake  boundary  con- 
ditions are  used  between  the  two  zones.   The  first  of  these 
demands  pressure  equilibrium  across  the  wake 


C^  =  C^8     at^so  (2.5.31) 


which  leads  to  the  equations 


fo^-jl^   =  ST  CC^  (2.5.32) 

w'3<«-Vi4(3Wirjjte>-y4re>]  =  FRUP      (2.5.33) 


where  the  right-hand  sides  of  (2.5.32,33)  refer  to  the 
steady  and  frequency  terms  respectively  of   C  ,   (2.3.9,10) 
The  second  wake  condition  requires  continuity  of  the  normal 
velocity  component  through  the  wake 


§£-  s  5£-      at«o  (2.5.34) 


which  leads  to  the  equations 


\ 


^CS)^^^)  so  (2.5.35) 


V»i(S^-*-h^-»-J±S!  ^iCS^V*1^]  =°       (2.5.36) 

Equations  (2.5.32,33,35,36)  are  solved  for  the  zone-Ill 
disturbances 

PL:  ^s(s^=0      h's(0)  =  -^-         (2.5.37) 

PI:  ^$^)=-jtS-0       VljlSisifS-^+U^       (2.5.38) 

2.5.5      Zone   V.       In   zone   V  an  unknown  disturbance   is   added 


(2.5.39) 


which  is  determined  by  enforcing  flow  tangency  at  the  lower 
surface  of  blade  two  in  zone  V 


PL:  9<(jt)-0  hsCo^rt-—-         (2.5.40) 


£s(.v)- O  hs(&)  =  --z- 


PI:  (2.5.41) 


which  leads  directly  to  the  pressure  difference  for  the 
last  portion  of  the  blade 


PL:        ACp^s-ZC^        »-m<^l        (2.5.42) 

PI:  (2.5.43) 

i-vn  <^  i  i 


2.5.6   Results  and  Conclusions.   The  finite  results  for 
the  second  blade  are  compared  with  four  existing  solutions 
in  figures  2.5.2  and  2.5.3  for  three  interblade  phase  angles: 
Kurosaka's  solution  [6]  is  a  closed-form,  analytical, 
infinite-cascade  solution;  Platzer's  [16,17]  is  a  numerical, 
method-of-characteristics  solution  for  the  finite  cascade; 
and  Adamczyk's  [19]  and  Verdon's  [12]  solutions  are  numeri- 
cal, infinite-cascade  results.   It  is  interesting  to  note 
that  the  finite  results  agree  exactly  with  Kurosaka's 
infinite  solution  over  the  entire  blade  for  the  plunge  case 
and  through  zone  II  for  the  pitch  case.   Over  the  zone-V 
portion  of  the  pitching  blade  —  the  portion  exposed  to  the 
wake  disturbances  from  blade  one  —  the  two  solutions  agree 
numerically,  as  seen  in  figure  2.5.3,  but  are  not  analytically 
the  same.   In  fact,  Kurosaka's  solution  is  singular  for 
M  =  0   (in-phase  oscillation  of  the  blades)  over  the  last 
portion  of  the  blade,  while  equation  (2.5.43)  behaves  well 
for  that  case.   All  five  solutions  were  found  to  be  in 
excellent  agreement  for  interblade  phase  angles  larger  than 
15°,  but  for  small  angles  considerable  disagreement  is 
evident,  as  seen  in  figure  2.5.2. 
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Figure  2.5.2  Comparison  of  the  finite-cascade 
solution  with  existing  solutions. 
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Chadwick  [5]  compared  his  finite  and  infinite  results 
for  the  pitch  case  over  the  in-flow  portion  of  the  blade  and 
noted  that  although  the  finite  and  infinite  pressure- 
difference  (AC  )   distributions  are  the  same,  the  local 
pressure  distributions   (C^,  CL)   may  in  fact  differ.   The 
author  found  this  phenomenon  occuring  over  the  remaining 
portion  of  the  pitching  blade  also.   Figure  2.5.4   shows  that 
the  local  pressure  distributions  for  the  present  finite 
solution  and  Verdon's  infinite  solution  do  not  agree  as 
well  as  the  pressure-difference  distributions.   Some  con- 
clusions can  be  drawn  from  the  analysis  at  this  point: 

1.  For  low  frequencies  of  oscillation  the  finite  cascade 
rapidly  approaches  the  infinite-cascade  solution,  achieving 
good  agreement  with  the  pressure-difference  distribution 

of  the  latter  at  the  second  blade. 

2.  For  small  interblade  phase  angles,  the  present  solution 
agrees  with  Kurosaka's  Laplace  transform  solution  for  the 
infinite  cascade  but  considerable  differences  with  the  other 
published  solutions  are  noted. 

3.  Although  the  present  finite-cascade  solution  and 
Verdon's  infinite-cascade  solution  for  the  pressure-difference 
distributions  are  in  agreement  for  large  interblade  phase 
angles,  the  local  pressure  distributions  resulting  from  the 
two  solutions  are  seen  to  disagree. 

2 . 6   The  Infinite  Two-Degree-of -Freedom  Cascade.   In  order 
to  demonstrate  the  applicability  of  the  above  analysis  to 
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Figure  2.5.4-  Comparison  of  the  finite  and  infinite 

local  pressure  distributions. 


the  infinite  cascade,  it  is  applied  to  an  infinite  version 
of  Cascade  B.   Since  the  advantage  of  an  identifiable  first 
blade  is  lost,  a  condition  specifying  the  space  periodicity 
of  the  flowfields  between  the  blades  is  derived.   The  pitch 
results  through  equation  (2.6.17)  were  first  obtained  by 
Chadwick  [5] ,    and  the  reflection  zones  referred  to  in  the 
solution  are  depicted  in  figure  2.6.1. 
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Figure  2.6.1   The  infinite-cascade  reflection  zones. 

2.6.1   The  Flow  Periodicity  Condition.   The  periodicity 
»f  the  flow  from  blade  passage  to  blade  passage  was  first 
.dentified  by  Kurosaka  [6]  and  is  the  heart  of  his  passage 
.pproach  to  the  infinite  cascade.   An  alternate  derivation 
ollows. 

Consider  the  corresponding  zones  I  and  IA,  and  assume 
.hat  the  potentials  are  the  same  there  except  for  a  complex 
actor  K 


9*«<*/^  *  K  SOS/O 


(2.6.1) 


iince  each  potential  must  satisfy  the  boundary  condition  of 
.ts  corresponding  blade 


iy-"^         mX^"°  (2.6.2) 


9  0H-W 


t 


-   -^e^  <x^  f  =0  (2.6.3) 


it   is   clear   that 


K-e~>  (2.6.4) 


and   so 


j6fcx,^  =  ^Cf^e-^  (2.6.5) 

The  same  result  follows  for  any  pair  of  corresponding  zones 
on  adjacent  blades  —  for  .pitch  and  plunge  —  and  is  termed 
the  flow  periodicity  condition. 

2.6.2.  Zones  I  and  IA.  The  analysis  in  zones  I  and  IA 
proceeds  much  as  in  the  finite  case  (section  2.5.1).  The 
velocity  potentials  have  the  form 

^  «£,<*)  +  *  [h,l&^^*^  *•<*:>]  (2.6.6) 

<?£>=&aM  +  k  [>W<rt-  ylJiA^  (2.6.7) 

and  the  boundary  conditions  yield 

PL:            ^«(S)-A     J)A(vr)-B>  (2.6.8) 


PI:  <Ji(s>*  J"^    c?^^>  =  ^  e^-hD      (2.6.9) 


hiCs-)  =  ^j  S  -  ^AS  +  E 
PL:  (2.6.10) 


The  constants   A   through   D   are  evaluated  using,  respectively, 
continuity  of  the  velocity  potential  across  the  leading-edge 
Mach  wave  of  the  (n+l)st  blade  and  the  flow  periodicity 
condition 

^(s^vn)  rr  9^(<3-  rO)  (2.6.12) 

0*«*S£  er^        (2.6.13) 

The  pressure  coefficients,  in  terms  of  those  of  the  isolated 
blade  (2.3.9,10),  follow 

PL:  CP°MX  =  C?,  (2.6.14) 


PI-  £$?    a  C  -  —  ■  B  •  x  (2.6.15) 


2.6.3   Zones  II  through  IV.   Finding  the  additional 
reflections  through  zone  IV  is  accomplished  by  routine  use 
of  the  flow  tangency  condition  (2.2.26,27),  and  so  only  the 
pressure  results  for  the  (n+l)st  blade  are  given  below 

PL:      ^Cp^'^x-iC^U-e^)  0<z^<Z^  (2.6.16) 

PI:  "  (2.6.17) 


O 


±  ^    <2.6<J 


PL:       ACp^-^  =  -^-C^C\-e^)  26<i<$<\-W       (2.6.18) 

PI:  ^  (2.6.19) 


2.6.4      Zones   V,    VI,    and  VII.      In  the    first   three    zones 
above   the  wake  of   the      nth  blade 

0*=.  $?+£s<s>+*  [hSW>-  vJ£^5c»J  (2.6.20) 

9&£  =  0?  *£♦<«+*  [h«OT-*  ^tj*«BJ  (2.6.21) 

*k  -  «,^w*^[M«-'f-v^3    (2-6-22) 

where  the   g.   and   h4   Sauer  functions  are  known  from  the 
analysis  in  zone  IV  and 


5s<*"~  V»sCO  =  o  (2.6.23) 

47C\-wo  =  V»7C»-vn)  =:  O  (2.6.24) 

2.6.5   Zones  VA  and  VIII.   In  the  zones  straddling  the 

wake  of  the  (n+l)st  blade,  using  flow  periodicity  and 
results  from  zones  IV  and  V 

*?«  -  tiZ  ♦,J|WkCrt+h[hw««- l^fl  ^Cr]  (2.6.25) 

-V-  1^5VisC0-)e";^+  W,C5=->  +  h»C^-)  (2.6.26) 


where 


(2.6.27) 


3*0)s  h%(0  =  O  (2.6.28) 


The  wake  boundary  conditions  -  pressure  equilibrium  and 
normal  velocity  continuity  -  are  enforced  between  zones 
VA  and  VIII 


°fl  m  C?*        ^    t-°  (2.6.29) 

■JY"   =   -JY"        **    t=°  (2.6.30) 


giving  rise   to   two   equations    in   three   unknown   Sauer 
functions 

jtU-€T*>  -£  V;  =  ST  [-*<($£-  O]  (2  .  6  .  31) 


^(t-e-^)+^;  +  /s  =  0  (2.6.32) 


(2.6.33) 


+  1 


[?s  o-  e^  -*  Vt]    =  fr  [-  itcC-  c{^ j| 


hVCwe^  +  w^^  +  ^^iO-e-^^+^+^e"^  =  O     (2.6.34) 

where  the  right-hand  sides  of  (2.6.31,33)  indicate  the 
steady  and  frequency  terms,  respectively,  of  the  enclosed 
quantity.   The  third  equation  in  the  set  comes  from  applying 
the  boundary  condition  to  <f>       ,   at  the  lower  surface  of 
the  (n+1) st  blade 


^e-^-^7  «  (2.6.35) 


V/5e-^-AT*  ^c^e-^-^-ysd/s  er*>  ^ 


=.0     (2.6.36) 


Solving  equations  (2.6.31)  through  (2.6.36)  for  the  Sauer 
functions  permits  determination  of  the  pressure  difference 
over  the  last  portion  of  the  blade 


PL: 


ACpe-^r-lC^      »-w<S  <2Lft+vn  (2.6.37) 


cPe-^  =  -ac£ 


-4-lfe   vw 


Pis  (2.6.38) 


l-vn<^<2.A«J-H'* 


PL:   ^Cpe-,•^=-2C;lC^-e:^)      2./}i+m<^±\  (2.6.39) 


PI:       -4-  [2y8«iMl-vC&^2KWi->.2^d)1o-e^)  -  /fr1^     (2.6-.40) 

2.6.6   Results  and  Conclusions.   Figure  2.6.2  shows  that, 
as  with  the  finite-cascade  solution,  there  is  considerable 
disagreement  among  the  several  existing  solutions  for  small 
interblade  phase  angles,  but  even  for  those  angles  the 
present  solution  and  Kurosaka's  are  in  agreement.   The 
pressure-difference  expressions  derived  above  for  the  infinite 
cascade  agree  exactly  with  those  of  Kurosaka  [6]  over  the 
entire  plunging  blade  and  up  to   f  st  Z.&A  ■*•  m.   for  the  pitching 
one.   Over  the  last  portion  of  the  pitching  blade  the  two 
solutions  are  slightly  different  analytically  but  are  in 
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Figure  2.6.2  Comparison  of  the  infinite-cascade 
'  solution  with  existing  solutions. 


KK 


numerical  agreement  as  seen  in  figures  2.6.3  and  2.6.4; 
the  analytical  discrepancy  could  not  be  resolved.   All  of 
the  solutions  are  in  accord  for  the  larger  phase  angles. 
Continuing  with  the  conclusions  begun  in  section  2.5.6,  it 
can  be  added  that: 

4.   For  interblade  phase  angles  larger  than  15°  the 
present  method  of  analysis  yields  closed-form  pressure- 
difference  expressions  which  are  in  agreement  with  existing 
solutions  for  the  slowly-oscillating  cascade. 

2. 7   Limiting  Cases  of  the  Finite  and  Infinite-Cascade 

Solutions . 

2.7.1   The  Steady  Cascade.   All  terms  of  the  plunge  solu- 
tions are  frequency  terms;  therefore  if   k   is  set  to  zero, 
the  pressure-difference  distribution  goes  to  zero  for  the 
entire  blade,  which  is,  of  course,  what  should  happen  for  a 
stationary  flat  plate  aligned  with  the  flow. 

If   k   is  set  to  zero  in  the  pitching  case,  then  one 
obtains  for  the  isolated  airfoil 


C?,^--|-  (2.7.1) 

r 


and  so  both  the  finite  (2.5.27)  and  infinite  (2.6.17) 
solutions  becomes 


=  __± (2.7.2) 
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Figure  2.6.3   Comparison  of  the  infinite-cascade 
solution  with  existing  solutions. 
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Figure  2.6.4   Comparison  of  the  infinite-cascade 
solution  with  existing  solutions. 


where  the  phase-angle  term  has  been  ignored  because  it  is 
meaningless  in  the  steady  case.   Equation  (2.7.2)  is  the 
classical  Ackeret  solution  [27]  for  a  flat  plate  at  an  angle 
of  attack  equal  to  one  radian,  which  is  the  maximum  up-pitch 
position. 

2.1.2  The  Single-Blade  "Cascade".   If  either  cascade  is 
"squeezed"  down  to  a  single  blade  by  setting  the  spacing 
parameters   m  and   d   to  zero,  then  the  solutions  (2.5.26,27) 
and  (2.6.16,17)  reduce  to  those  of  the  isolated  blade 

ACP=  -2.C£  (2.7.3) 

where  the  phase  term  has  again  been  ignored. 

2.7.3  The  Infinitely-Spaced  Cascade.   The  isolated-blade 
solution  is  again  recovered  from  equations  (2.5.27)  and 
(2.6.17)  as  the  cascade  spacing  is  allowed  to  get  large. 

If  the  spacing  parameters   m  and   d   approach  infinity  in 
such  a  way  that  the  terms  ySdM2"   and   ftzm  remain  equal, 
then  equation  (2.7.3)  again  results.   The  significance  of 
these  two  terms  is  discussed  further  in  the  next  section. 

2.8   The  Unique  Inflow  Condition.   A  consideration  of  the 
infinitely-spaced  cascade  in  section  2.7.3  leads  one   to 


The  author  is  grateful  to  Professor  M.  F.  Platzer  for 
first  asking  this  question. 


<<:   For  which  cascades  does  the  relationship 


*<jM*=/a*m  (2.7.4) 


ad?   As  seen  from  equation  (2.5.27),  for  such  a  cascade 
:e  pressure-difference  over  the  inflow  portion  of  the  blade 
i  merely  that  of  the  isolated  pitching  blade  times  the 

fictor  CI-  e"^). 

Let  o<  be  the  Mach  angle  defined  by 


u^lN-'(^  (2«7'5) 


hen  it  is  not  difficult  to  show  that  for  any  cascade  having 
he  stagger  angle 

,  ...  /  sin 2.*   \         (2.7.6) 

!quation  (2.7.4)  will  also  hold.   Said  another  way,  all 
„enfcers  of  an  infinite  family  of  pitching  cascades  with  the 
stagger  angle  given  by  (2.7.6)  -  and  with  no  restrictions 

,  .n  h-,ve  exactly  the  same  inflow  pressure 
on  blade  spacing  -  will  have  exactly 

distribution,  i.e. 


PI: 


^Cp  =  -*c?,u-e--)  (2-7.7) 


This  is  termed,  by  the  author,  the  unique  inflow  condition 
for  the  slowly-oscillating  cascade.   Figure  2.8.1  depicts 
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Figure  2.8.1   The  unique  inflow  condition 


three  different  cascades,  each  with  the  same  stagger  angle, 
satisfying  (2.7.6).   The  pressure  distribution  for  each 
was  obtained  via  Platzer's  method-of-characteristics  pro- 
gram.  The  unique  inflow  condition  is  evidenced  by  the 
exact  agreement  of  the  pressures  for  all  three  cascades 
over  the  first  portion  of  the  blade. 

2.9   Some  Comments  on  Resonance.   For  a  given  reference 
blade  in  an  infinite  cascade,  resonance  occurs  when  aero- 
dynamic disturbance  waves  generated  by  upstream  blades  are 
in  phase  with  the  motion  of  the  reference  blade  [12] . 
Samoylovich  [2  8]  shows  that  at  the  resonance  condition  the 
linear  cascade  experiences  zero  aerodynamic  damping.   The 
development  that  follows  is  based  on  Verdon ' s  [12]  excellent 
explanation  of  the  resonance  phenomenon  and  provides  a  means 
of  delineating  what  he  calls  the  subresonant  and  superresonant 
regions  of  the  u-te  plane. 

Referring  to  figure  2.9.1,  also  taken  from  [12],  it  is 
seen  that  the  leading  edge  of  each  blade  creates  disturbances 
which  propagate  in  all  directions  at  the  freestream  speed 
of  sound   a^ .   These  disturbances  are  swept  downstream  by 
the  freestream  flow  U  ,  and  there  are  exactly  two  resultant 
waves  which  propagate  along  the  cascade  leading  edge  in 
the  direction  opposite  to  the  direction  of  rotation  and  at 
velocities   V,   and  V?    .   From  the  law  of  cosines  one  has 
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Figure  2.9.1   The  backward-traveling  waves 
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Figure  2.9.2   The  subresonant  and  superresonant  regions 


The  times   tln#  t2n   required  for  a  wave  from  an  upstream 
blade  (n  =  -1,  -2,  .  ..)  to  reach  any  given  reference  blade 
(n  =  0)  are 

^n/in=^[M^^)±[vn(W^2^J/2|        (2.9.2) 

The  reference  blade  lags  an  upstream  blade  by 


"to-  "»ff-       Wa-»ra;...  (2.9.3) 


and  resonance  occurs  when  the  upstream  waves  are  in  phase 
with  the  reference-blade  motion,  i.e.  when 

tosU,»  .  (2.9:4) 

or  when 

/VsZ~p"  ^tVn-»-Bci)±[m(w-v.2.ya<i)],/a}       (2.9.5) 

For  a  fixed  reduced  frequency   k   the  superresonant  inter- 
blade  phase  angles  lie  between  mk     and   M2 /  as  shown  in 
figure  2.9.2,  while  the  subresonant  phase  angles  lie  to 
either  side.   The  superresonant  region  is  important  because 
Verdon  is  of  the  opinion  that,  in  that  region,  neither  a 
finite-cascade  solution  nor  a  passage-approach  infinite- 
cascade  solution  will  adequately  approximate  the  true 
infinite-cascade  solution.   In  an  effort  to  shed  some  light 


on  the  subject,  some  comparisons  are  presented  in  figures 
2.9.3  through  2.9.6.   For  cascade  A,  with   k  =  0.05  ,  the 
superresonant  range  of  phase  angles  is  approximately   -7° 
to   -2°  ,  and  so   -5°   was  examined.   Figures  2.9.3  and  2.9.4 
show  that  Platzer's  finite  solution  and  the  author's  passage- 
approach  infinite  solution  are  in  agreement  with  Verdon ' s 
solution.   For  cascade  B,  however,  Verdon 's  solution  is 
significantly  different  from  the  other  two  over  the  last 
portion  of  the  blade  for  the  superresonant  phase  angle   -3°  . 
Figures  2.9.5  and  2.9.6  show  some  high-frequency  comparisons 
for  the  phase  angle   -90°  ,  which  is  in  the  superresonant 
regions  of  both  cascades  A  and  B.   In  these  two  cases  the 
two  solutions  —  one  finite  and  one  infinite  —  are  in  accord. 
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Figure  2.9.3   Low-frequency  superresonant 
pressure  distributions. 
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pressure  distributions. 


3.   THE  NONLINEAR  CASCADE 

3.1   General.   In  two  papers  Teipel  provides  the  numerical 
tools  for  the  analysis  of  low  supersonic  flow  past  a  single 
oscillating  airfoil  of  finite  thickness.   In  one  paper  [23] 
he  presents  an  unsteady,  nonlinear  characteristics  method 
to  determine  the  continuous  variation  of  flow  variables. 
In  the  other  [24]  he  treats  the  oscillating  shock  waves 
present  in  the  nonlinear  problem  to  permit  determination 
of  the  discontinuous  jumps  in  the  flow  variables  across  the 
shocks  emanating  from  the  single  blade.   Using  these  tech- 
niques as  well  as  some  computational  and  theoretical  insight 
provided  by.  Chadwick  [5] ,  who  very  ably  translated  and  applied 
Teipel' s  work,  the  present  study  is  aimed  at  finding  the 
pressures,  steady  and  unsteady,  on  the  second  blade  of  a 
parabolic-arc  cascade.   The  cascade  oscillates  in  pitch 
alone,  and  the  pressures  are  sought  over  the  entire  upper 
blade  surface  and  over  the  lower  surface  up  to  the  point  of 
impingement  of  the  trailing-edge  shock  from  the  preceding 
blade  (figure  3.2.2).   These  surface  pressures  are  computed 
by  means  of  a  FORTRAN  IV  computer  program  prepared  by  the 
author  for  the  IBM- 360  computer.   The  program  listing  appears 
at  the  end  of  this  work  immediately  after  Appendix  A.   What 
follows  is  an  orderly  presentation  of  the  theory  underlying 
the  computer  program  and  the  computational  techniques  used 
in  implementing  that  theory.   Wherever  possible  the  text 
notation  will  coincide  with  that  used  in  the  program. 


3. 2   Cascade  Geometry. 

3-2.1   Blade  Shape.   The  blade  shape  was  chosen  to  approx- 
imate the  blades  in  the  test  cascade  employed  by  Fleeter 
and  Riff el  [25]  who  attempted  to  measure  steady  and  unsteady 
pressures  in  a  pitching  cascade.   The  airfoil  used  in  the 
computer  model  has  a  flat  upper  surface  and  a  lower  surface 
conforming  to  the  parabolic  arc  given  by 


df 


=  4-t  *  fX-\)  (3.2.1) 


where  chordwise  distance   x  and  midchord  thickness  X.     are 
normalized  by  chordlength.   Figure  3.2.1  shows  the  computer 
model  for  t  =  0.03   superimposed  on  Fleeter' s  test  blade 
and  lists  the  surface  coordinates  of  the  test  blade  in  per- 
cent chord.   The  greatest  discrepancy  between  the  two  shapes 
occurs  on  the  aft  portion  of  the  blade  in  the  region  behind 
the  lower-surface  shock  reflection;  however,  the  lower- 
surface  pressures  are  not  computed  in  this  region. 

3.2.2   Blade  Spacing.   The  blade  spacing  parameters  are 
input  variables  for  the  program  and  are  discussed  in  more 
detail  in  section  3.5.1.   The  physical  spacing  of  the  blades 
of  the  computer-model  cascade  —  dubbed  Cascade  T  —  is  shown 
in  figure  3.2.2. 

3.3   Teipel's  Method  of  Characteristics.   The  underlying 
insight  of  Teipel's  treatment  of  the  unsteady  problem  is 
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Figure  3.2.1   Test  blade  shape  and  coordinates 
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Figure  3.2.2   Geometry  of  the  computer  model. 

i  representation  of  the  velocity  potential  as  the  sum 
I  a  steady  component  and  a  small -amplitude  -unsteady  one 


flfatftfcl  =  &c*/2>  +  A<*/?>^L 


kt. 


(3.3.1) 


:  depicted  in  figure  3.3.1.   When  (3.3.1)-  is  substituted 
Jto  the  transonic  unsteady  small-disturbance  equation 


=  0   (3.3.2) 


id  into  the  unsteady  boundary  condition 


1For  a  discussion  of  the  limitations  imposed  by  this 
Juation  see  Teipel  [23]  and  Chadwick  [5]. 


Figure  3.3.1 
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Figure  3.3.2   The  characteristic  grid. 


S£   ~  SX   ^  at  (3.3.3) 


where   y(x,t)   describes  the  oscillating  airfoil  surface, 
one  obtains  two  separate  problems  coupled  by  the  steady 
velocity  potential  0      .   One  problem  is  nonlinear,  steady 


and  zero-angle-of-attack 


and  the  other  linear  and  unsteady 


[/**«*&]3l "  7?>  p^lfe + «^-mWA- 


(3.3.4) 


(3.3.5) 


O       (3.3.6) 


Sfc-sa  -  Tl-*-  Lte<X-b)l  (3.3.7) 

9^;  L  J 


where   0  (x,y)   is  an  amplitude  function  and  the  factor  ec 
is  omitted.   It  is  important  to  note  that  since  the  char- 
acteristic directions  are  determined  by  the  coefficients  of 
the  highest  derivatives,  both  problems  will  have  the  same 
characteristic  directions.   The  state  variables   for  the  two 


1The  symbol  /x  is  used  both  for  the  steady  state  variable 
and  the  interblade  phase  angle,  and  it  is  normally  clear 
from  usage  which  meaning  applies. 


problems  are  essentially  the  steady  and  unsteady  velocity 


components   u  ,  v  ,  u..   and  v. 


STEADY:  ^^+MZ(UI)U0  (3.3.8) 


y*:=3/2M20r  +  0V*  (3.3.9) 


UNSTEADY:  U|,  V»  (3.3.10) 


where 


fia=M2-l  (3.3.11) 

The  characteristic  directions  for  both  problems  are • 

^-^]        =  ■+"  -i-  (3  3  12) 

where  the  upper  sign  refers  to  the  left-running  ©<  charac- 
teristics and  the  lower  to  the  right-running  A     character- 
istics.   The  compatibility  relations,  however,  are  not 
the  same  for  the  two  problems.   They  are 

STEADY:          7^/z  X  M=  CONST^  4  (3.3.13) 


A^X  yUr:  COtslST^ 


1  2       -*- 

The  symbol  yfi  is  used  both  for  the  quantity    (M  -l)2 
and  the  right-running  characteristic  direction.   It  is 
normally  clear  from  usage  which  meaning  applies. 


UNSTEADY:  (3.3.14) 


where 


^«s|^-^cl*     ApsJX  (3.3.15) 


Equation  (3.3.14)  requires  some  additional  comment. 
Note  first  that  in  addition  to  the  state  variables   u, 
and   v,  ,  0,   is  also  present  in  the  unsteady  compatibility 
relations.   Since 

«**■?§■*' - -if  <•;>    '  (3-3-16' 

then   0,   must  be  found  by  integration  as  one  moves  from 
one  gridpoint  to  another.   For  the  points  A,  B,   and   C 
shown  in  figure  3.3.2,  for  example 

«:  0,c=0*  +  \(u»-».  -\)<A*  (3.3.17) 

ys:  $*f  =  j6*  +  j«*i--^jW*         (3.3.18) 

In  performing  the  actual  computations  an  average  of  the  two 
above  expressions  is  employed. 


The  factor  ^  in  (3.3.14)  is  straightforward,  repre- 
senting the  increment  in  *   along  the  B      characteristic. 
At  point  C  in  figure  3.3.2  for  example  one  could  use  the 
approximation 


» 


^/  =  *c-*fc  (3.3.19) 


However  z^^    is  a  bit  more  complex,  for  while  the  increment 
dx   is  along  the  <x  characteristic,  the  quantity  $2±S\ 
represents  the  change  in  X  with  respect  to  x  along  the 
crossing  B     characteristic  at  the  point  of  interest. 
Again  at  point  C  one  could  employ 


^^  =  ^~a  CXc-*^  (3.3-.20) 


3.4   Teipel's  Treatment  of  the  Oscillating  Shock  Wave  with 
Extension  to  the  Present  Cascade. 1 

3.4.1   The  Isentropic  Rankine-Hugoniot  Pressure  Relation. 

Consider  a  normal  shock  moving  through  an  undisturbed 

perfect  gas  at  velocity  V   relative  to  the  gas  as  shown 

in  figure  3.4.1,  where  the  hat  notation  indicates  conditions 

behind  the  shock.   The  governing  continuity,  momentum  and 

energy  equations  are 


Except  for   u,,  v.,  lu   and   v", ,  the  velocities  appearing 
in  sections  3.4.1  through  3.4.6  are  dimensional,  in  a  brief 
departure  from  the  normalization  of  all  velocities  by  the 
freestream  velocity   U  . 
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Figure  3.4.1   Jump  conditions  across  a  moving, 

normal  shock. 
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Figure  3.4.2   Consequences  of  the  isentropic 

assumption. 


^pC-Vs^  =r^tV*-V&-V^  (3.4.1) 


P+y>Vs2  =  P  -v^CV^-Vs-V^  (3.4.2) 


h+^  s  ^^CV^Vs-V^  (3>4#3) 


Assuming  constant  specific  heats,  the  speeds  of  sound  are 
introduced,  and  the  enthalpies  are  eliminated  by 


a.z=-*i-      0^=*-^  (3.4.4) 


h= K  =  3.4.5) 


and,  after  a  good  deal  of   manipulation,  the  change  in 
the  normal  velocities  across  the  shock  is  found  in  terms 
of  the  shock  velocity  and  sound  speed  ahead  of  the  shock 


%*-■&-<*--*>  (3.4-6) 


Equations  (3.4.1)  and  (3.4.6)  and  the  isentropic  relation 


-M-^r 


yield  the  Rankine-Hugoniot  pressure  relation  sought 


p     y-n[u'  J 


(3.4.8) 


Figure  3.4.2  is  based  on  figure  8-3  of  [29],  and  it  shows 
that  the  isentropic  approximation  introduces  no  significant 
error  for  Fleeter 's  test  cascade  or  its  computer  model. 

3.4.2   The  Blade-Two  Shock  Position.   The  leading-edge 
shock  from  the  pitching  blade  oscillates  out  of  phase  with 
the  blade.   Let  tf  (y)   be  the  shock  angle  of  the  steady, 
leading-edge  shock  as  shown  in  figure  3.4.3.   Since  the 
shock  is  curved,   y   varies  with  distance   y   from  the 
blade,  and  by  convention 

0<^0<^-  (3.4.9) 


The  shock  angle  at  any  point  on  the  oscillating  wave  can 
be  written  as  a  small  variation  from  the  steady  one 


V  =s  *0  +X'elk*  (3.4.10) 

In  order  to  locate  the  oscillating  shock  relative  to  the 
steady  one,  Teipel  introduces  the  complex  displacement 
function   g(y)  ,  also  shown  in  figure  3.4.3,  and  so  the   x 
coordinate  of  a  shock  point  is  approximately 

■X-  'Xo  ■+  (>-<^  COT*o  +  Jty  e."^      (3.4.11) 


1The  symbol  gamma  used  here  and  in  figures  3.4.4  and 
3.4.5  represents  the  total  shock  angle.   All  other  uses  of 
the  unsubscripted,  unsuperscripted  gamma  stand  for  the 
dimensionless  ratio  of  specific  heats. 
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Figure  3.4.3  The  oscillating  leading-edge  shock  wave, 
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where   (X0'Y0)   are  the  steady  coordinates  of  the  blade's 
leading  edge.   Taking  the  total  differential  of  (3.4.11) 
and  noting  that  along  the  shock 


=  -  COT  tf  (3.4.12) 


=># 


results  in 


*  v  i^_  cwW-t 


S  =  &IN'^0^-  fe*-***  (3.4.13) 


thus  relating  shock  angle  and  shock  displacement. 

3.4.3   The  Velocity  Components  Upstream  of  the  Blade-Two 
Shock  Wave.   Ahead  of  the  shock  the  velocity  components 
have  both  steady  and  unsteady  parts 

w.*  =  U  +  ufl  +  u,  e.1**  (3.4.14) 

V*  s  v0  -v  v,  e'c**  (3.4.15) 

where  the  barred  velocities  are  dimensional  counterparts  of 
the  dimensionless  quantities   u  ,  v  ,  u.^   and  v^   already 
introduced  in  equations  (3.3.8,9,10).   However,  the  flat  upper 
surfaces  of  the  blades  of  cascade  T  are  aligned  with  the  flow 
in  the  steady  position,  and  so,  upstream  of  the  shock 

U0=Vo-=:0  (3.4.16) 


Referring  to  figure  3.4.4  one  has  for  the  normal  and 
tangential  velocity  components  ahead  of  the  shock 

Vn-=TJ  SIN  tf0  +  U  COS  Vo  tt'  -*•   (  U,  SlNtfo  +  V»  COS  *o}  ^""^  (3.4. 17) 

Vt^Ucos^o-XI  s\n*»*'  -*•  (H,  cos*0~  \7,  sin  *m  felv,,t    (3.4  .18) 

3.4.4   The  Blade-Two  Shock-Motion- Induced  Velocity. 
Relative  to  the  gas  ahead  of  it  the  shock  is  moving  at 
the  velocity 

Vs^-Vk+v'  (3.4.19) 

where   V   is  a  normal  velocity  induced  by  the  motion  of 
the  shock,  and  is  seen  from  figure  3.4.5  to  be 


V'^l^-  2>IN*  (3.4.20) 


or,  from  (3.4.11) 

V'  =  U  ^iNtfo  9  e>Wt  (3.4.21) 

3.4.5   The  Velocity  Jump  Conditions  Across  the  Blade-Two 
Shock.   Using  equation  (3.4.13)  to  eliminate   tf'   from  the 
expressions  for  the  upstream  velocity  components  (3.4.17,18), 
they  become 
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Fiaure    3.4.5      The   shock-motion-induced  velocity. 


Vtts"U  sin  Vo  f  (U  costfo  sm^o^  u»  &iN*.+vlcos*o)e'-kt'     (3  . 4  .  22) 

"a 


V^Ucostfo-K-'UsiN^oi^-  ■*■  u,cos«6-v,  ^N^)elw   (3.4.23) 


Note  that  each  of  the  above  velocities  has  three  component 
parts.   The  first  term  in  each  expression  arises  from  the 
steady  shock  position,  the  second  arises  from  the  oscillatory 
motion  of  the  shock  about  that  steady  position,  and  the 
third  and  fourth  terms  arise  from  the  unsteady  disturbances 
ahead  of  the  shock  which  are  generated  by  the  first  blade's 
motion. 


The  components  behind  the  shock  are  found  by  noting 
that  the  tangential  component  does  not  change 


V^=V^  (3.4.24) 


and  by  using  (3.4.6,19,21,22)  to  obtain 


<U^LU    ^0O+^^) 


+  Rfer^  cos*°     *iNl*o  U-TTT7I>4^-  (3.4.25) 


where 

MN=M    S\NVo  (3.4.26) 

It  is  convenient  to  once  again  return  to  the  dimensionless 
x-y  velocity  components  u,  and  V,  ,  obtaining  from  the 
unsteady  terms  of  (3.4.24,25) 


2 


Ui  -=.  wi,-r^-  -*■  <■  W*  9  -+  miu>-*-  vn^v,        (3.4.27) 


%  •=.  Y\x-r-    "*■  ^i^  +  ^u,  -m^v,         (3.4.28) 
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where   the   factor      e      c     has   been  dropped   and 


m'=--^T    SVN   **•    SlNl^0  (3.4.29) 


m**  ^^T  KCi"*"mS'>    *»Nz*0  (3.4.30) 


vnj=cosz^  -v  ~^-    ^n**,^--^*.  )  (3.4.31) 


^= ~  +  VTT— —  U-7T   mT)  (3.4.32) 


K»  =  —  ~-(ws2>f8+  -—-)    &in*V%  (3.4.33) 


Ha*  -^—  kCl-*-  -V^    COTtfo    SIN***  (3.4.34) 


V^Ytt^.  (3.4.35) 


K4.  =  sin**q  +  ^-7  co«aV«  U--~ Jr-)  (3.4.36) 


Note  that  the  above  coefficients  are  functions  of  the 
local  shock  angle  "tf    and  so  are  not  constant  along  the 
shock . 

If  no  unsteady  disturbances  exist  upstream  of  the  shock 
(u,  =  v.  =0)   then  (3.4.27,28)  reduce  to  Teipel's  single- 
blade  results.   If,  in  addition,  the  shock  does  not  oscillate 
(g  =  dg/dy  =0)  ,  then   u,   and  v,   are  identically  zero 
indicating  that  there  are  no  unsteady  disturbances  behind 
the  shock.   Equations  (3.4.27,28)  represent  an  extension  of 
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Teipel's  work  from  the  single  blade  to  the  cascade,  and 
they  were  first  obtained  by  Chadwick  for  a  cascade  of 
wedges.   It  is  noted  for  completeness  that  the  above 
equations  would  apply  equally  well  to  an  upper-surface 
shock  wave  if  the  signs  of   m,  ,  m,  ,  n_   and  n,.   were 
reversed. 

3.4.6   The  Pressure  Jump  Conditions  Across  the  Blade-Two 
Shock.   An  additional  extension  of  Teipel's  work  was  accom- 
plished by  the  author  in  order  to  develop  a  more  accurate 
method  of  determining  the  flow-variable  discontinuities 
than  that   employed  by  either  Chadwick  or  Teipel.   The 
improved  method  is  treated  in  more  detail  in  section  3.8, 
but  it  is  based  on  the  extension  of  the  Rankine-Hugoniot 
pressure  relation  (3.4.8)  to  the  cascade  via  -equations 
(3.4.19,20,21),  resulting  in 


?-ar'^-  +  L^+f»tt,+  M     (3-4-37) 


where 

J^=  ^~-    M*      S1N£X0  (3.4.38) 

k    =--±^-  m£  fe  (3.4.39) 

J?3=£t^n  (3-4-40' 
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*+9TZT    M^  C°T  *°  (3.4.41) 


Equation  (3.4.37)  relates  the  pressure  jump  across  the 
shock  to  the  shock  displacement  and  the  upstream  unsteady 
velocities.   In  the  single-blade  case   u,   and   v,   are 
zero,  and  Teipel's  result  is  recovered. 

3.4.7   The  Boundary  Conditions  for  the  Shock  Equations. 
In  section  3.8,  the  flow  variables  behind  the  shock  are 
obtained  by  the  simultaneous  solution  of  the  difference 
equations  resulting  from  (3.3.14)  and  the  shock  equations 

(3.4.27,28,37).   In  order  to  start  the  solution  process, 
then,  several  quantities  must  initially  be  known  at  the 
point  of  intersection  of  the  shock  and  blade,  i.e.  at 

(x  ,y  ).   From  figure  3.4.3  it  is  clear  that 

o<u0->  ^O  (3.4.42) 

Also,  immediately  behind  the  shock  the  flow  tangency 
condition  for  the  second  blade  must  hold 

0»s-[l-«-  LkCX0-b>]e^   A-fc  fXo,y^      (3.4.43) 

With  v,   known  equation  (3.4.28)  provides 

S^- J-^-*  Ul-V^V,)     <k±    CXo»       (3.4.44) 
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and  (3.4.44)  used  with  (3.4.27,37)  gives  the  remaining 
boundary  conditions 

*»S*TV»  + Wx u>  +     n>   *  V>    *fc  llU,^.)     (3.4.45) 


fo©  Too  hi  rli  Vli 


(3.4.46) 


a-t  LKo^o) 


Of  course  for  equations  (3.4.43)  through  (3.4.46)  the  hatted 
quantities  are  just  behind  the  shock  at   (x  ,y  )  . 
Equation  (3.4.45)  was  derived  by  Chadwick  while  (3.4.46) 
is  the  author's  result. 

3.4.8   The  Velocity  Potential  Jump  Conditions  Across  the 
Blade-Two  Shock.   Landahl  [26]  shows  that  the  total  velocity 
potential  (3.3.1)  is  continuous  across  the  leading-edge 
shock,  but  this  does  not  imply  that  either  the  steady  or 
unsteady  component  is  individually  continuous.   Teipel  and 
Chadwick  choose  to  assume  that,  at  the  first  shock  point 


#iCX©,#©)^  ^.CXc^o)  (3.4.47) 


The  author,  on  the  other  hand,  chooses  to  find  the  initial 
velocity  potential  jump  at  the  point   (x0»yo^  *   For  a 

perfect  gas  with  constant  specific  heats,  it  can  be  shown  that 
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ft--v<**tk  (3'4-48> 

and  since  u\  and  P/P»  are  known  at   (xQ,yo)   from 
equations  (3.4.45,46),  then  so  is  \    .   Figure  3.4.6  is  a 
plot  of  computer  data  for  cascade  T,  and  it  shows  the 
discontinuity  in  ^  along  the  shook  as  a  function  of 
distance  from  the  blade.   Note  that  the  initial  discon- 
tinuity is  small  justifying  Teipel's  assumption  (3.4.47), 
but  equation  (3.4.48)  eliminates  having  to  make  that 
assumption  at  all. 

3  5  Constructing  the  Characteristic  Grid.   Because  the 
characteristic  directions  are  the  same  for  both  the  steady 
and  unsteady  problem,  the  characteristic  grid  can  be  con- 
structed and  employed  to  find  both  the  steady  and  unsteady 
velocities.   The  grids  above  blade  one  (zone  A,  and  above 
blade  two  (zone  C)  are  composed  of  two  families  of  straight, 
parallel  characteristics.   The  mesh  below  blade  two  (zone  B, , 

however,  is  composed  of  one  straight,  divergent  family   (/> 

and  one  curved  family   («0  • 

3  51  The  Straight  Character!^  -  Grids  Above  B_lades_One 

^^o.   since  the  upper  blade  surfaces  generate  no  steady 

i-x    -x    q\    +-hat,  in  zones  A  and  C 
1  +.  .;  G  copn  from  (o.o.o)    un^1-/  ->-" 
disturbances,  it  is  seen  "^ 

KA*l-fi*  (3.5.1) 
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Figure  3.4.6   The  discontinuous  jump  in  the  unsteady 
velocity  potential  across  the  leading- 
edge  shock. 
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and  so  the  characteristics  in  those  zones  are  simply  Mach 
waves  (figure  3.5.1).   Each  blade  has  21  equispaced  upper- 
surface  gridpoints,  and  the  program  requires  that  the  leading 
edge  of  the  second  blade  be  positioned  at  one  of  the  straight- 
mesh  gridpoints.   This  is  accomplished  via  the  input  parameters 
K   and   N  which  fix  the  cascade  blade  spacing  according  to 


YB  =  O.Oi5  T  (3.5.2) 

r 


X&  =  Y3  fi  ■+■  O.OS"  N  (3.5.3) 

For  cascade  T,  K   and  N   are   16   and   7   respectively. 
Figure  3.5.1  also  shows  the  characteristic  numbering  conven- 
tion employed  in  the  program,  so  that  a  gridpoint  may  be 
located  by  two  characteristic  coordinates 


a,;^)  <3-5-4> 


Note  that  the  straight  mesh  is  computed  as  if  both  blades 
were  flat  plates.   The  x-y  coordinates  of  the  straight 
mesh  are  referred  to  as   XI,   Yl   in  the  computer  program. 

3.5.2   The  Curved  Characteristic  Grid  Below  Blade  Two. 
In  zone  B  (figure  3.5.2)  the  characteristic  grid  is  non- 
linear and  must  be  constructed  iteratively.   Figure  3.5.2 
reveals  the  numbering  convention  employed  for  the  curved 
mesh.   The  curved-mesh   x-y   coordinates  are  referred  to  as 
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Figure  3.5.2  The  curved  characteristic  grid. 
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X2,Y2   in  the  program  to  distinguish  them  from  those  of 
the  straight  mesh.   Before  proceeding  with  the  grid  con- 
sturction,  some  useful  relations  are  obtained. 

For  an  <X  characteristic  passing  through  the  shock 
from  zone  A  to  zone  B  (figure  3.5.3) 


>A  -/A  ^  >%-/**  (3.5.5) 


and  since  u*    is  zero 

*a=  t*A*+/*Oy3  (3.5.6) 

Now  consider  two  points  on  the  same  jS     characteristic, 
also  illustrated,  in  figure  3.5.3.   From  the  <X  compatibility- 
relation 

?v(X,3"/2-  aa{-L,X)  =  7^  (3.5.7) 

*lX,VO%-MCX,K)=***a  (3.5.8) 

and  from  the  A   compatibility  relation 

MX,CD*2-*.  u (!,<!)  a  MX,K)%  -v-^CX^K)  (3.5.9) 

and  so  it  must  be  that 

>CX/T)=MX,K)     A4<X^)=/*tt,yO  (3.5.10) 
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which  proves,  in  addition,  that  the  A      family  is  composed 
of  straight  lines. 

The  steady  shock  angle  at  any  point  along  the  shock  is 
found  by  averaging  the  angles  of  the  characteristic  lines 
before  and  after  the  shock,  i.e. 

*°=T  [taN-,(X-*WTAN-|(V4)]        (3.5.11) 

Taking  the  tangent  of  (3.5.11)  and  applying  the  sum  and 
double-angle  trigonometric  formulae  leads  to  a  quadratic 
equation  in  tan  tf 

C,XV*+/X,4)TAN2-30  +  ZCKV*^V*-I)TAN  a*  -  (*''*+ ^v*)=0     (3.5.12) 

The  solution  of  (3.5.12)  is  rather  cumbersome  and  is 
adquately  approximated  by  the  simpler  expression 


TAN  **:=   -,Z  >,,.  (3.5.13) 


In  the  computational  implementation  of  (3.5.13),  referring 
to  figure  3.5.3  again,  the  characteristic  angles  at  adjacent 
points  behind  the  shock  are  first  averaged  together,  and 
the  result  is  averaged  with  the  zone-A  angle,  i.e. 

TAN*«:r  — »  ;  — ^rr   —   .,  (3.5.14) 


2£V*  +  Ml^'i+M^K)''* 


_£7_ 


3.5.2.1  Computing  the  Point  (1,1).   The   X2-Y2   coord- 
inates of   (1,1)   are  of  course  those  of  the  blade-two 
leading  edge.   /A (1,1)   is  found  from  the  steady  boundary 
condition  (3.3.5)  and  equation  (3.3.9)  while  7\(1,1)   follows 
from  (3.5.6) . 

3.5.2.2  Computing  the  Points  (1,2)  and  (2,2).   The 
points   (1,2)   and   (2,2)   are  established  iteratively  as 
follows.   Program  subroutines  employed  are  indicated  in 
parentheses,  and  figure  3.5.4  is  provided  for  reference. 

1.  Establish  the  reference  point  XA  on  the  blade, 
close  to  the  leading  edge 

XA.=  XZ(M)  +  o.oos 

2.  Compute  u    at   XA  ,  which  is  equal  to  U(l,2) , 
using  the  steady  boundary  condition  (3.3.5)  and  (3.3.9). 

3.  Compute   A  at   XA  ,  which  is  equal  to  7^(1,2), 
using  equation  (3.5.6). 

4.  (Subroutine  SOLVE)  Locate   (1,2)   at  the  intersection 
of  a  line  through   (XA,  Y2(l,l))  with  slope 


TUi,^* 


and  a  line  through   (X2(l,l),  Y2(1,D)  with  slope 

2± 

2/?^  +^»/i),/l  +   M»,ttv* 
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Computing  curved-mesh  points  (1,2)  and  (2,2). 
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shock  wave 
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<XAfY2(l,l)')  ^2,2) 
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Computing  subsequent  shock- front  points  (1,1) 


Fiqure  3.5.4   Computing  the  shock-front  points. 
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5.  (Subroutine  ITER)  Guess  at  the  position  of   (2,2) 
on  the  blade 

7J.tztz^  •=   xzcy) 

6.  (Subroutine  ITER)  Increment  that  guess 

XT(2,2)  =  X2(2,2.)  -t-COOl 

7.  (Subroutine  ITER)  Compute  JU{2,2)       from  (3.3.5,9). 

8.  (Subroutine  ITER)  Compute  A(2,2)   from  (3.5.6). 

9.  (Subroutine  ITER)  Solve  for   X2(2,2)   from  the 
slope  of  the  <X  characteristic  through   (1,2)   and   (2,2) 

X2U,2)-XZO,Z)  Z 


10.  (Subroutine  ITER)  Does   X2(2,2)  =  X2(2,2)? 

YES:   GO  TO  11;  NO:   GO  TO  6 . 

11.  STOP. 

3.5.2.3   Computing  Subsequent  Shock  Front  Points. 
Subsequent  points   (1,1)   along  the  shock  are  also  found 
iteratively.   The  technique  follows,  and  figure  3.5.4 
applies. 

1.   Relocate   XA  on  the  blade 


where   DXA   is  based  on  the  first  position  of 


position  of   XA 


DXA=  XI(2,Z)-XA 


2.   (Subroutine  SHOKPT)  Guess   A(1,I) 


3.  (Subroutine  SHOKPT)  Increment  that  guess. 

Ml,X)  =  MlyX^  +  O.OOI 

4.  (Subroutine  SHOKPT)  Compute  U(1,I)      using  AA   and 
A  at   (2,1)  ^ '  which  are  the  same  respectively  as  M(2,2) 
and  7v(2,2) 

5.  (Subroutine  SOLVE)  Locate   (1,1)   at  the  intersection 
of  a  line  through   (XA,  Y2(l,l))  with  slope 


JxL 


Ml;XV/2" 


and  a  line  through   (X2 (1,1-1),  Y2 (1,1-1))   with  slope 

-4- 
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6.  (Subroutine  SHOKPT)  Compute  yu(l,l)   from  the 
steady  boundary  condition  (3.3.5)  at   XA   and  equation 
(3.3.9) . 

7.  (Subroutine  SHOKPT)  Does  U{1,I)    =  Z<(1,I)? 
YES:  GO  TO  8;  NO:  GO  TO  3. 

8.  STOP. 

3.5.2.4   Computing  the  Curved-Mesh  Interior  Points.   A 
point   (I, J)   interior  to  the  curved  mesh  is  located  (figure 
3.5.5)  at  the  intersection  of  a  line  passing  through 
(I-1,J)  with  slope 


and  a  line  through   (1,1)   with  slope 


-1 


3.5.2.5   Computing  Subsequent  Blade-Surface  Points. 
Subsequent  points  of  the  form   (1,1)   on  the  lower  surface 
are  computed  in  a  manner  exactly  analagous  to  the  computation 
of  the  point   (2,2)   as  described  in  steps  5  through  11  of 
section  3.5.2.2. 

3.6   Computing  the  Steady  Velocities  and  Pressures.   With 
U   and  7<    known  throughout  the  grid,  the  steady  velocities 
follow  directly  from  (3.3.8,9).   The  steady  pressures  on 
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Figure  3.5.5   Computing  the  curved-mesh 
interior  points . 


the  upper  surface  of  blade  two  are  zero  and  on  the  lower 
surface  are  given  at  the  surface  points   (1,1)   by 


Co  -  -  Z.  Uo  cx/x^ 


(3.6.1) 


3.7   Computing  the  Unsteady  Flow  Variables  at  the 
Straight-Mesh  Grid  Points. 

3.7.1   The  Blade-One  Upper-Surface  Grid  Points.   In  the 
straight  mesh  the  unsteady  velocities  and  velocity  potential 
are  designated  by 


ux  =  UllR  +  i  UllI     vL  =  VllR  +  i  VllI 


fa   =  FlR  +  i  FlI 


(3.7.1) 


The  unsteady  boundary  condition  for  blade  one  (3.3.7) 
allows  determination  of  the   v,   velocities 


VllR( 


1,1)  =  -1   V11I(I,I)  =  -k[xl(I,I)  -  b]     (3.7.2) 


Using  the  j$      compatibility  relation  (3.3.14)  in  difference 
form  results  in  a  system  of  linear  equations  for  the 
unknowns   U11R(I,I),   0111(1,1)  .   Figure  3.7.1  A  applies, 
and  the  equations  appear  in  matrix  form  in  section  A.l  of 
the  appendix.   The  velocity  potential  is  found  by  integration 
of  the  velocities,  using  equation  (3.3.18);  see  section  A. 2 
of  appendix  A. 

3.7.2  The  Straight-Mesh  Field  Points.   As  seen  in  figure 
3.7.1  B,  grid  points  not  on  the  blade  require'  use  of  both 
compatibility  relations.   The  equations  resulting  for  the 
unknown  u..   and  v,   velocities  at  each  field  point   (I/J) 
are  given  in  A. 3  of  the  appendix.   Equations  (3.3.17,18) 
are  averaged  to  find  the  velocity  potential;  see  section 

A. 4  of  the  appendix. 

3.7.3  The  Grid  Points  Immediately  Upstream  of  the  Shock. 
Figure  3.7.1  C  shows  that  the  shock  points  lie  in  the  straight 
mesh  and  therefore  the  flow  variables  at  each  shock  point  — 
upstream  of  the  shock  —  are  found  by  means  of  an  interpolation 
algorithm.   The  distances   RI   in  figure  3.7.1  C  are  computed 
and  any  desired  quantity   Q(1,K)   is  found  from  the  values 

of   Q   at  the  four  surrounding  points 
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Figure  3.7.1   Computing  the  straight-mesh  and 
curved-mesh  grid  points. 
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Bn  (-3.  7.  5. 


RD 


where 


4- 
RRTT^yfRX        Rfc=lll  RR.7        (3.7.4) 


3 . 8   Determining  the  Jumps  in  the  Flow  Variables  Across  the 
Blade-Two  Shock.   It  is  in  the  determination  of  the  jump 
conditions  that  the  author  employs  a  technique  —  alluded  to 
in  section  3.4.8  —  that  is  different  from  the  method  used 
by  Teipel  in  the  single-blade  case.   The  essence  of  Teipel's 
approach  is  to:  • 

1.  Assume  that  the  velocity  potential  does  not  change 
across  the  shock  at  the  first  point  (3.4.47). 

2.  Use  the  shock  velocity  equations  (3.4.27,28),  with 
boundary  conditions  (3.4.42,43,45),  in  conjunction  with  the 

&    compatibility  relation  (3.3.14),  to  solve  for  u1#  v± 
(behind  the  shock)  and   g   at  each  shock  point. 

3.  Use  equation  (3.3.18)  behind  the  shock  to  recover 
0        at  each  shock  point. 

In  contrast,  the  method  employed  in  the  present  study  is  to: 

1.   Use  (3.4.46,48)  to  determine   0,   behind  the  shock 
at  the  first  shock  point. 


106 


2.  Use  the  shock  velocity  equations  (3.4.27,28),  the 
shock  pressure  equation  (3.4.37),  the  &     compatibility 
relation  (3.3.14)  and  the  boundary  conditions  (3.4.42,43, 
45,46)  to  solve  for   i^,  'v^,  g   and   P/p   at  each  shock  point. 

3.  Use  equation  (3.4.48)  to  find  0        from  the  pressure 
jump  at  each  shock  point. 

The  latter  method  results  in  a  larger  system  of  equations 
at  each  point,  but  no  assumptions  concerning  the  continuity 
of  the  velocity  potential  are  required. 

3.8.1   Details  of  the  Jump-Condition  Computations. 
Step    2   above  requires  some  additional  comment.   The  shock 
equations  (3.4.27,28,37)  are  applied  in  difference  form  from 
one  shock  point  to  the  next.   They  are  the  last  six  equations 
in  the  matrix  equation  given  in  appendix  section  A. 5.   The 

£   relation  supplies  the  first  two  equations  in  that  matrix 
equation.   Note  from  figure  3.7.1  D  that  using  the  £     relation 
at  shock  point  (1,1)  requires  two  tacit  assumptions.   First, 
the  flow  variables  at  the  point  where  the  &      characteristic 
through   (1,1)   leaves  the  blade  are  essentially  those  at 
the  point   (1,1)   immediately  behind  the  shock.   Second,  the 
flow  variables  on  the  Q     characteristic  through   (1,1+1) 
and  just  behind   (1,1)   are  about  the  same  as  those  at 

(1,1)   behind  the  shock. 

In  the  computer  program  the  flow  variables  at  shock 
point   (1,1)  /  ahead  of  the  shock,  have  the  form 

uL  =  U12(1,I)  -  i  U12I(1,I)  (3.8.1) 


while,  at  the  same  point  immediately  behind  the  shock, 
they  are  designated 

uL  =  URJ(I)  +  i  UIJ(I)  (3.8.2) 

until  stored  in  the  form  of  (3.8.1). 

3. 9   Computing  the  Unsteady  Flow  Variables  at  the  Curved- 
Mesh  Grid  Points. 

3.9.1  The  Blade-Two  Lower-Surface  Grid  Points.   The 
unsteady  boundary  condition  for  the  second  blade  provides 
the  lower-surface   v,   velocities  (figure  3.9.1  A).   The 
equations  are  given  in  section  A. 6  of  the  appendix.   The 

u,   velocities  are  found  via  the  <=<    compatibility  relation 
(section  A. 7)  and  the  velocity  potential  follows  from 
equation  (3.3.17)  (appendix  section  A. 8). 

3.9.2  The  Curved-Mesh  Field  Points.  Both  compatibility 
relations  are  used  at  a  field  point,  as  seen  from  figure 
3.9.1  B.   The  resultant  matrix  equation  appears  in  appendix 
section  A. 9.   The  velocity  potential  is  found  from  equations 
(3.3.17,18)  (appendix  section  A. 10). 

With  the  curved  mesh  completely  computed  the  unsteady 
pressures  at  blade-two,  lower-surface  grid  points  are  given 
by 

Cn,  =  -2[ul2R(I,I)  -  k  F2I(I,I)| 
Pl      r  i.  (3.9.1) 

-2i[ul2I(I,I)  +  k  F2R(I,I)J 
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Figure  3.9.1   Computing  the  curved-mesh  grid  points. 

3.10   Results  and  Conclusions. 

3.10.1   The  Isolated-Blade  Case.   If  the  unsteady  boundary 
condition  along  the  first  blade  is  set  to  zero,  the  second 
blade  sees  freestream  conditions  and  its  unsteady  pressures 
can  be  compared  with  Teipel's  isolated-blade  results  [23]. 
This  is.  accomplished  in  figure  3.10.1,  where  it  is  noted 
that  a  small  discrepancy  between  the  two  solutions  exists. 
This  discrepancy  is  attributable  in  part  to  the  different 
techniques  used  by  Teipel  and  the  author  in  solving  the 
shock  equations  (see  section  3.8).   But  it  is  also  noted 
that  Chadwick  [5]  too  encountered  a  discrepancy  when  he 
applied  Teipel's  method  to  the  single,  biconvex  airfoil. 


109 


Figure 


3.10.1   is 


olated-Blade  Results 


♦■that  for  the  single  blade,  as  k 
Teipel  points  out  that,  ror 

_*,  >-   This  phenomenon  is  also  ev 

blade  in  cascade  «.  .^  results  of 

PXatzer's  linear  solution  [171 

_       v   of   1 • 5  • 
this  investigation  for  a 

„  i  in  3  shows 
«,.  SteadyPressure^esults,   Figure  3.10.3 

3.10-2  The_Stead2_ —  two_ 

for  the  lower  surtace 
the  steady  pressures  experimental  results 

It  ig  seen  that  agreement  between  the  exp 


110 


-4.0 


Cascade 

T 

M: 

1.55 

blade  2 

lower  surface 

k:  1.5 

A 

30.0° 

Platzer- 

-P 

Strada-S 

0.1 

> 


0.2 

1 


0.3 
1 


0.4 


Figure  3.10.2   High-k  comparison  of  the  linear 

and  non-linear  solutions. 


of  Fleeter  [25]  and  the  numerical  results  of  the  computer 
model  is  good  at  midchord  and  fair  at  the  leading  edge. 
The  program  results  compare  favorably  with  theoretical 
predictions  from  second-order  Prandtl-Busemann  theory  and 
a  third-order  theory  due  to  Ferri  [30]  which  takes  entropy 
losses  into  account.   This  latter  comparison  again  demon- 
strates that  the  computer  solution  is  not  limited  by  its 
use  of  the  isentropic  Rankine-Hugoniot  pressure  relation 
(3.4.8) . 

3.10.3   The  Unsteady  Pressure  Results.   Figures  3.10.4 
through  3.10.12  show  the  computer  results  of  the  present 
study,  the  experimental  data  from  Fleeter 's  investigation, 
and  the  linear  results  of  Platzer's  method-of -characteristics 
program  for  nine  interblade  phase  angles.   The  pressure 
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Figure  3.10.3   Blade-two  lower-surface 
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data  are  presented  as  an  amplitude  and  phase  angle  relative 
to  the  motion  of  the  second  blade.   The  phase  angle  is 
measured  from  the  point  where  the  second  blade  attains 
maximum  up  pitch  to  the  point  of  peak  blade  pressure.   Along 
the  upper  surface  of  blade  two  the  linear  and  nonlinear 
numerical  results  agree  exactly  -  because  the  upper  surface 
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is  flat  -  providing  an  additional  check  of  the  nonlinear 
program.   Along  the  lower  surface  the  nonlinear  effects 
show  up  as  differences  between  the  two  predicted  pressure 
distributions  —  differences  which  vary  with  interblade  phase 
angle.   The  predicted  nonlinear  pressure  amplitude  is,  in 
all  cases  examined,  at  least  as  large  as  the  linear  predic- 
tions, which  is  as  expected.   In  figure  3.10.5  (m=  -4.4°) 
the  linear  and  nonlinear  amplitude  predictions  are  nearly 
coincident,  while  in  figure  3.10.9  (M=  130.6°)  they  differ 
considerably.   Figure  3.10.13  is  a  plot  of  the  difference 
between  the  nonlinear  and  linear  amplitude  predictions  versus 
the  turning  angles  experienced  by  the  flow  at  the  leading 
edge  of  blade  two.   It  shows  that  the  nonlinear  computer 
predictions  .a"re  consistent  since  the  variation  from  the 
linear  solution  increases  with  the  flow  turning  angle. 
That  turning  angle  is  given  by 

T.Av.=  57.  3°  U-  COS^W)  -t-  6.8°  (3.10.1) 

when  blade  two  is  in  the  maximum  up-pitch  position.   Note 
that  the  large  turning  angles  in  figure  3.10.13  result  from 
the  choice  of  a  one  radian  maximum  pitch  amplitude  for 
nondimensionalization  purposes. 

Comparison  of  the  nonlinear  computer  predictions  with 
Fleeter 's  experimental  data,  however,  reveals  that  the 
theory  and  the  experimental  results  do  not  consistently  agree 
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Figure  3.10.4   Blade-two  unsteady  pressures 
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Figure  3.10.5   Blade-two  unsteady  pressures 
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Figure  3.10.6   Blade-two  unsteady  pressures 
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Figure  3.10.7   Blade-two  unsteady  pressures 
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Figure  3.10.8   Blade-two  unsteady  pressures 
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Figure  3.10.9   Blade-two  unsteady  pressures 
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Figure  3.10.10   Blade-two  unsteady  pressures 
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Figure  3.10.11   Blade-two  unsteady  pressures 
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Figure  3.10.13   Nonlinear  variations  from  linear  theory. 

In  figure  3.10.11,  for  example,  there  is  adequate  upper- 
surface  agreement  with  the  theory,  but  the  lower-surface 
pressure-amplitude  comparison  is  poor.   In  fact,  in  no  case 
is  there  agreement  between  data  and  theory  for  both  phase 
and  amplitude  along  the  lower  surface.   The  lower-surface 
experimental  data  appear  to  be  somewhat  inconsistent.   They 
sometimes  vary  widely  from  one  point  to  the  next  for  the 
same  interblade  phase  angle  (figure  3.10.6)  and  they  also 
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vary  a  great  deal  for  rather  small  changes  in  interblade 
phase  angle  (figures  3.10.10,11,12).   m  addition,  in  each 
of  the  figures  there  is  at  least  one  amplitude  data  point 
for  the  nonlinear  (lower)  surface  that  is  below  the  corres- 
ponding linear-surface  data  point.   This  phenomenon  is 
contrary  to  expectation  because  the  thickness  surface  must 
generate  a  larger  magnitude  disturbance  than  the  flat  surface. 
Certainly,  additional  experimental  data  are  required  to 
adequately  evaluate  the  theoretical  methods. 

3.10.4   Conclusions.   In  summary,  it  is  concluded  that  the 
extension  of  Teipel's  single-blade  methods  to  the  nonlinear 
cascade  results  in: 

1.  Steady  pressure  predictions  which  are  in  accord 
with  other  nonlinear  theories  and  which  are  supported,  to 
some  extent,  by  experimental  results. 

2.  Unsteady  pressure  predictions  which  exhibit  expected 
trends  relative  to  linear  theory,  but  which  are  not  consis- 
tently in  agreement  with  experimental  data  presently  available. 
The  significant  influence  of  thickness  on  the  unsteady 
pressures  is  apparent  in  both  the  theory  and  data. 

3.  A  numerical  approach  to  the  unsteady  cascade  which 
appears  to  have  the  potential  to  treat  the  effects  of  blade 
thickness,  blade  camber,  and  steady  angle-of -attack  loading 
to  provide  a  more  realistic  model  of  the  fan  and  compressor 
blades  in  a  turbomachine. 
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APPENDIX  A 


A.l  Equations  for  the  Unknown  u-\    Velocity  at  Blade-One 
Surface  Points, 

1.  References:  section  3.7.1;  figure  3.7.1  A. 

2.  Auxiliary  constants: 

AX    =   2kM2(0.025)/7vA 
A3    =    1+  0.0125   A2 
A5   =   1-0.0125  A2 

3.  Matrix  equation: 


A2    =  k2M2(0.025)/?vA 
A4    =   A3/7^ 
A6    =   A5/7y£ 


A5  -Al/2 

Al/2  A5_ 


U11R(I,I) 
U11I(I,I) 


A3   U11R(I-1,I)  + 
A3    U11I( 1-1,1)  - 


♦  Al    UlII(I-l,I)/2  -   A4   VllR(I,I)+A6   VllR(I-l,I)+ 

-_A1    UllR(I-l,I)/2  -   A4   VllI(I,I)+A6   VllI(I-l ,  I)  + 

+  A2    F1R(I-1,I) 
+  A2    F1I(I-1,I) 


A02   Equations    for  &    at   Blade-One   Surface   Points. 

1.  Reference:    section  3.7.1;    figure   3.7.1   A. 

2.  Equations: 

F1R(I,I)    =    FlR(I-l,I)    +{ullR(I,I)   +  U11R(I-1,I)    - 

-    [V11R(I,D  ^  V11R(I-1,I)]^}(0-0125) 

F1I(I,I)    =    F1I<I-1,I)    +{ullI(I,I)    +  13111(1-1,1)    - 

-   [yiii(i,i)  +  viii(i-i,i)]/^/lV0-0125) 


125 


A.  3  Equations  for  the  Velocities  un  and  v-|  at  Straight -Mesh 
Grid  Points. 

1.  References:  section  3.7.2;  figure  3.7.1  B. 

2.  Auxiliary  constants: 

Al    =    2kM2(0.025)/>,  A2    =  k2M2(0.025)/^ 


A3    r   1   -   0.0125  A2 
A5    =   A3/xV* 

3.  Matrix  equation: 


A4   =   1  +  0.0125   A2 
A6   =   A4/>Yi 


A3 

-Al/2 

-A6 

Al/2 

A3 

0 

A3 

-Al/2 

A6 

Al/2 

A3 

0 

0 

■A6 
0 
A6 


U11R(I,J) 
U11I(I,J) 
VllR(I,J) 
V11I(I,J) 


A4  U11R(I,J-1)  + 
A4  UllI(I,J-l)  - 
A4  U11R(I-1,J)  + 
A4    U11I(I-1,J)   ~ 


+  A1   U11KI,  J-l)  -  A5   V11R(I,J-1)  +  A2    FlR(I,J-l) 

Z 

-  Al    U11R(I,J-1)  -  A5  V11I(I,J-1)  -V  A2    FlI(I,J-l) 

Z 

+  Al   UllI(I-l,J)  +  A5  V11R(I-1,J)  -V  A2    FlR(I-l,J) 
2L 

-  Al.  U11R(I-1,J)  -*-  A5   V11KI-1,  J)  +  A2    FlI(I-l,  J) 


A. 4    Equations   for  jr\_   at    Straight -Mesh   Grid  Points. 

FlR(I,J)    =  i[FlR(I,J-l)  +  FlR(I-l,J)]    +   0.0125   UllR(I,J)    ■+ 
+  0.00625  |ullR(  I,  J-l)   +  UllR(I-l,J)]    -*■ 
+  0.00625  [yllR(I, J-l)    -   VllRCl-ljJ)]/^4 

FlI(I,J)    =   i[FlI(I,J-D  -+  Fl  1(1-1,  J)]     +  0.0125    U11I(I,J)    + 
+  0.00625  [ylll( I,  J-l)    +  UllI(I-l,J)]     -f 
+  0.00625  [vllKl,  J-l)    -   V11I(I-1,J)]  A* 
1.    References:    section  3.7.2;    figure   307.1   B. 
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A. 5  Equations  for  the  Jump  Conditions  across  the  Blade-Two 
Leading-Edge  Shock. 

1.  References:  section  3.8.1;  figure  3.7.1  D. 

2.  Auxiliary  constants: 

Al    =   2/  LMl,I)  +7U1,I-1)J 

A2   =   4kM2[x2(l,I)  -   X2(1,I-1)]/  [>d,D  +7^(1,1-1)] 

A3    =   A2  k/2 

DY  =     [Y2(1,I)    -   Y2(l,I-l)]/2 

3.  Matrix  equations: 


1 

-pa/2 

Al 

0 

-A3/2 

0 

0 

0 

IRJ(I) 

A2/2 

l 

0 

Al 

0 

-A3/2 

0 

0 

UIJ(I) 

1      • 

0 

0 

0 

0 

0 

-m1/DY 

Tn~ 

VRJ(I) 

0 

l 

0 

0 

0 

0 

-m^ 

-m1/DY 

VIJ(I) 

0 

0 

1 

0 

0 

0 

-nL/DY 

n2 

FRJ(I) 

0 

0 

0 

1 

0 

0 

-n2 

-nx/DY 

FIJ(I) 

-*N2 

0 

0 

0 

0 

w\ 

-Pi/DY 

P2 

GR(D 

0 

-1M2 

0 

0 

-*M2k 

0 

-P2 

-PL/DY 

GI(I) 

URJ(M)  +  A2  UIJCI-D/2  ♦  Al  VRJ(I-l)  +  A3  FRJ(I-l)/2  • 
UIJ(I-l)  -  A2  URJ(I-D/2  +  Al  VIJ(I-l)  -t  A3  FIJ(I-l)/2 
m3[ul2R(l,I)  +  U12R(1, 1-1)]  *  ™4  [vl2R(l,I)  +  Vl2R(l,I-l)] 
m3[ul2I(l,I)  +  U12I(1,I-1)]  ^m4[vl2T(l,I)  +  Vl2I(l,  1-1)] 
m4[ul2R(l,I)  -*U12R(1,I-1)]  *  n4  [vl2R(l,  I)  +  Vl2R(l,  1-1)] 
m4[ul2I(l,I)  +  U12I(1,T-1)]  +  n4  [vl2I(l,  I)  +  Vl2I(l,  1-1)] 
<P3  "  M2)[U12R(1,I)  ♦  U12R(1,I-1)]  +  p4[vl2R(l,I)  + 
|(P3  -  M2)[U12I(1,I)  +  U12I(1,I-1)]  +  p4[vl2I(l,I)  + 
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A. 5  (Continued^)  : 

3.  Matrix  equation  (right  hand  side  continued): 


-  URJ(I-l)    -  m1    GR(I-1)/DY  -  m2   GI(I-l) 

-  UIJ(I-l)    -  m-L   GKI-D/DY  -vm2  GR(I-l) 

-  VRJ(I-l)    -  n1   GR(I-1)/dY  -   n2   GI(I-l) 

-  VIJ(I-l)    -  nL   GKI-D/DY  +  n2   GR(I-l) 

♦  V12R(I-1)]    +  k*H2[F2I(I)  +  F2I(I-1  )]  +X  M2URJ(I-1 )  - 
+  V12KI-1)]    -  k*M2[F2R(I)  +  F2R(I-1)]**M2UIJ(I-1)  + 


-*M2k   FTJ(I-l)    -   p1   GR(I-1)/DY-  p2   GI(I-l) 
-»-*M2lc   FRJ(I-l)    -   px    GI(I-1)/DY  +  p2    GR(I-l) 


A. 6   Boundary-Condition  Equations    for  the  v-|    Velocity  at   the 
Blade-Two   Surface   Points. 

1.  References:    section  3.9.1;    figure  3.9.1  A. 

2.  Auxiliary  constants:      X  =  X2(I,I)    -  X2(l,l) 

3.  Equations: 

V12R(I,I)    =    -    COS(yU)    -  k(X-b)SIN(yu) 
V12I(I,I)    =    -    SINOO   -   k(X-b)COS(ytf) 
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A. 7  Equations  for  the  ui  Velocity  at  Blade-Two  Surface 
Points, 

1.  References:  section  3.9.1;  figure  3.9.1  A. 

2.  Auxiliary  constants  (also  used  in  A. 9  equations): 
DXA  =  X2(I,J)  -  X2(I-1,J)     XLA  =  i[>(I,J)  +  X(I-l,'jJ] 
DXB  =  X2(I,J)  -  X2(I,J-1)    XLB  =  |[X(I,J)  +  A(I,J-1)] 
VAR  r  -|-[vl2R(I, J)  +  Vl2R(I-l,J)J   DLB  =  Xl,J)  -  X(I,J-1) 
VAI  =  i[vl2I(I,J)  +   V12I(I,J-1)]       Al  =  2  k  M2 

DVAR  =  Vl2R(I,J)  -  V12R(I-1,J)        A2  =  k  Al/2 

DVAI  =  Vl2I(I,J)  -  V12I(I-1,J) 

SLA  =  i  DV(I,J)V  XCI-lfJ)^  SLB  =  i£>(I,j/%  X(I,J-l!P!] 

DLXB  =  [>(I,J+1)  -  Ml,J)]/[x2(I,J+l)  -  X2(I,J)] 

3.  Matrix  equation: 

1 


DLXB  DXA  .  A2  DXA2     _  Al  DXA 
4  XLA      2  XLA 

Al  DXA 


2  XLA 


1  ♦ 


2  XLA 

DLXB  DXA   A2  DXA2 
4  XLA      2  XLA 


U12R(I,J) 
U12I(I,J) 


U12R 
U12I(I 


1(1-1, J)[l 
-l,J)[l 


DLXB  DXA  4.  A2  DXA21  ,  ltl,T/T  i  T^  Al  DXA 
4  XLA   +  2  XLA  J  +  Ul2I(I  1>J)   2  XLA  * 

A2^2]-  Ul2R(I-l,J)  A1_DXA^ 
2  XLA  J  '2  XLA 


1  -  DLX3  DXA 
4  XLA 


DVAR    A2  VAR  DXA2 


SLA 


XLA 


SX3-  «*"  F2R(I-1,J)  -3 


A2  DXA 


XLA 


DVAI    A2  VAI  DXA2  A2  DXA 

SLA  "*   XLA  SLA   «+  F2I(I-1  ,  J)  XLA 


4.  Note:  at  surface  points  I  =  J. 
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A. 8  Equations  for  fa    at  Blade-Two  Surface  Points. 

1.  References:  section  3.9.1;  figure  3.9.1  A. 

2.  Auxiliary  constants:  as  defined  in  A.7. 

3.  Equations: 

F2R(I,J)  r  F2R(I-1,J)  -V- ji[ul2R(I,  J)  *  Ul2R(I-l,  J)] 

+  var/sla^  dxa 

F2I(I,J)  =  F2I(I-1,J)  +   £|[ul2I(I,J)  +  Ul2I(I-l,J)] 

•v  vai/sla}  DXA 

4.  Note:  at  surface  points  I  =  J. 

A. 9  Equations  for  the  Velocities  ut  and  v^  at  Curved -Mesh 
Grid  Points . 

1.  References:  section  3.9.2;   figure  3.9.1  B. 

2.  Auxiliary  constants:  as  defined  in  A. 7; 

3.  Matrix  equation: 


I    .  DLXB  DXA  _  A2  DXA2    _Al  DXA   „  1  ,.  A2  DXA2       0 
4  XLA      2  XLA       2  XLA    SLA   2  XLA  SLA 


Al  DXA   ,  . .  DLXB  DXA  _  A2  DXA2     0     1  _,  A2  DXA2 

2  XLA     mr    4  :CLA      2  XLA  SLA   2  XLA  SLA 


DLB  _  A2  DXB2    _A1  DXB     1  .l.  A2  DXB2         0 
4  XLB     2  XLB        2  XLB    SLB    2  XLB  SLB 


Al  DXB    ,     DLB   .  A2  DXB2        n  1   ,  A2  DXB2 


1     DLB  __  A2  DXB^        n        1 
i"4  XLB  ~   2  XLB         U      "SLB-*" 


2  XLB     ^  4  XLB     2  XLB  SLB  ^  2  ](LB  SLB 

(Note:  continued  on  next  page.) 
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A. 9    (Continued): 

3.   Matrix  equation   (unknown  vector  and  right -hand- side 
matrix  continued): 


U12R(I,J) 


U12I(I,J) 


Vl2R(I,J) 


V12I(I,J) 


U12R(I 


U12I(I 


-1    j)    [~i        DLXT 
L,J;    L  4   > 


-It J)    [}   ' 


DLXB    DXA    .    A2 
XLA  2 


DLXB    DXA  ^  A2 


DXA2] 

xla  J 


4  xla 


U12R(I, J-l) 


U12I(I,J-1) 


DLB  .  A2    DXB2 


4   XLB 


DLB 
4   XLB 


DXA2"|      _ 
XLA  J 


] 


2   XLB 


A2 


DXB2]      _ 
XLB   J 


^,a2ia.itj)JLm^vi2aa.ifj)  [-^  +  J£^]  + 

-ni2R(i.i,j)4^^vi2i(i-i,j)  [-^+j«j#£[|* 

+  U12I(I,J-1)4^,V12R(I,J-1)  [^SLb'2AXLBXSLb]^ 

.UI2R(I,J.1)^^.*V12I(I,J-1)  [^-i^fft]^ 


F2R(I-1,J)    ^    PXA 
XLA 


+  F2I(I-1,J)^DXA 
XLA 


+  F2R(I,J-D    ^    PXB 


+  F2I(I,J-1)    ^y8 
XLB 
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A. 10  Equations  for  fa    at  Curved-Mesh  Grid  Points, 
L.  References:  section  3.9.2;  figure  3.9.1  B. 

2.  Auxiliary  constants:  as  defined  in  A. 7. 

3.  Equations: 

F2R(I,J)    =   -§-[f2R(I-1,J)  +  F2R(I,J-1)] 

+  i[ul2R(I,J)  +  Ul2R(I-l,J)  ^V12R(I7J)  +  Vl2R(I-l,  J)]pXA 

+  i{lJ12R(I,J)  +  U12R(I,  J-l)   -  V12RCI.J)    -  V12R(I,  J-l)lnYP 
*  <.  '  SLB  -> 

F2I(I,J)    =  i[F2I(I-l,J)  +  F2I<I,J-1)] 

+  i?Ul2I(I,J)  +  U12KI-1, J)  ^  V12I(I,J)  +  yl2I(I-l,J)WA 

*>  SLA  J 

*  i(ul2I(I,J)  +  Ul2I(I,J-l)   -  V12I(I7J)^V12I(I,J-1)] 


DXB 
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